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Abstract 

The equations of classical spin-orbit motion can be extended to a 
Hamiltonian system in 9-dimensional phase space by introducing 
a coupled spin-orbit Poisson bracket (jjlj) and a Hamiltonian 
function (|20|), After this extension and by establishing connections 
between initial and extended systems it becomes possible to apply the 
methods of the theory of Hamiltonian systems to the study of 
polarized particle beam dynamics in circular accelerators and storage 
rings. Some of those methods have been implemented in the com- 
puter code FORGET-ME-NOT |, g. 
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1 Introduction 



We began writing this paper in the summer of 1994 year following the 
suggestion of D.Barber to summarize our results of the investigation of the 
dynamics of unpolarized and polarized proton beams, which we obtained 
during 1991-1993 and which were only partly available and only in the very 
compressed form of conference proceedings. 

After starting this work it became clear that it would be better to sep- 
arate the available material into several papers, and for the first one we 
chose those results which can be incorporated together with help of the re- 
search technique which we call Hamiltonian extension of the equations 
of classical spin-orbit motion. Such important topics as the computa- 
tion and analysis of the one-turn Taylor maps for spin and orbit motion, 
spin dependent invariant functions and their connection with the Derbenev- 
Kondratenko vector, the normal form algorithm for equations of spin motion 
in the SU(2) representation and, of course, many more practical subjects 
related to the question how to preserve the polarization of a particle beam 
during acceleration will not appear in the present paper and will be published 
later. 

For various reasons the publication has been prolonged for almost three 
years, and now we are glad to be free from this, but, at the same time, 
we are not too enthusiastic about the necessity to write up the remaining 
unpublished results. 



1.1 The Classical Equations of Spin-Orbit Motion 



The quasi-classical description of the motion of a relativistic nonradiat- 
ing point particle with spin in accelerators and storage rings includes the 
equations of orbit motion which we write in the Hamiltonian form 

dq _ dH orbt dp _ _ dH orb t 

dt dv 



and the Thomas-BMT equation 

ds 
dt 

Here 



dt dq 
for the classical spin vector s 



(1) 



W x s 



(2) 



H 



orbt 
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G ( 7T • B ) 7T 

W = — I (1+jG)B 2 V 9/1 , — I G 



mojc \ mQC 2 (l + 7) m c y 1 + 7, 



n x £ 



and £ is the time. The vectors q = (91,92,93) and p = (^1,^2,^3) 
are canonical position and momentum variables, and s — (si, S2, S3) is the 
classical spin vector of length h/2. The parameters e and mo are the charge 
and the rest mass of the particle, c is the velocity of light, G — (g — 2)/2 
which quantifies the anomalous spin g factor, 7 is the Lorentz factor, 7? 
is kinetic momentum vector, £ and £> are the electric and magnetic fields, 
and A and $ are the vector and scalar potentials. 

B = curl 9 A (3) 
1 cM 

£ = -grad g -$ - - — (4) 



c dt 



e r 

7T = p — - A 

C 



H orbt - e$ 
T = 2 — 



\ 



1 + p- (5) 
V m c/ 



Later on we will refer to the system (TO-(|2J) as the triangular system (the 
equations of spin motion contain the orbital variables but the evolution of 
orbital variables does not depend on the spin degree of freedom). 

Remark: In the Hamiltonian picture of orbital motion we cannot define 
the Lorentz factor (§), in the usual way, in terms of the absolute value v of 
the particle velocity v 

1 2 . . 

7 = , V = V ■ V 

\A - (v/cf 

but instead need to express it through canonical variables. 
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2 The Hamiltonian Extension of the Equa- 
tions of Classical Spin-Orbit Motion 



The Poisson bracket lies at the basis of the modern point of view of the 
Hamiltonian formalism. The idea of axiomatic introduction of the bracket 
very likely belongs to Dirac. In this section we recall the necessary definitions 
from the theory of Hamiltonian systems and introduce coupled spin-orbit 
Poisson brackets which give us the canonical extension of the equations of 
classical spin-orbit motion. 

2.1 Poisson Bracket 

Let M be a finite-dimensional manifold (phase space) and C°°(M) be 
the linear space of smooth functions f : M —> R . The binary operation 

{*, *} : C°°(M) x C°°(M) -> C°°(M) 

called the Poisson bracket satisfies the relations 

bilinearity : {Xf + fig, h} = X {/, h} + fi {g, h} 

antisymmetry : {/, g} = -{g, /} 

Leibnitz rule : {/, g ■ h} = {/, g} ■ h + g ■ {/, h} 

Jacobi identity : {{/, g}, h} + {{h, /}, g} + {{g, h}, /} = 

where A, /x are arbitrary constants. 

Let £ be local coordinates on a manifold M . From the bilinearity and 
Leibnitz rule it follows that for any fixed function h , the linear operator 



is a differentiation and hence may be represented in the form (see, for example 



{h, *} : C°°(M) 



C°°(M) 



1,0) 



{h, *} = At(z) 



d * 

dzi 



(6) 
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(here and further on the summation over repeated indices is implied). From 
(H) one finds by direct substitution of functions Z{ that 

A h t (z) = {h, Zi } (7) 

Comparing (|GJ) and (0) we obtain 

{h, *} = {h, Zl } |^ (8) 

Consider the Poisson bracket {/, g} . By successively applying (||) to the 
expressions {/, *} and {*, Zj} we obtain 

Now we introduce the skew-symmetric matrix J = ({z^, zj}) and represent 
the Poisson bracket @ in the form 

{/, 9} = grad^/- J grad^p (10) 

Thus in a fixed local coordinate system z , the Poisson bracket is completely 
defined if we know the values = {z iy Zj} as functions of z . 

The matrix J is called the structural matrix of a Poisson bracket. 



2.2 Hamiltonian Dynamical Systems 

Again let z be local coordinates on a manifold M . Hamiltonian 
systems by definition have the form: 

^ = & H ^ ( U ) 

where H = H{r, z) is an arbitrary function (possibly depending explicitly 
on tq) , called the Hamiltonian. Using the representation of Poisson bracket 
([T0| ) we can rewrite equations (JTTJ) 



as 



dz 

- = J(z)gmd s H (12) 



1 Until now the equations of classical spin-orbit motion (0)-(||) are written using the 
time t as independent variable, but later on we will change this independent variable to 
be a path length along the design orbit z. This is the reason for us, here and further on, 
to formulate general results, which are not sensitive to the specific form of the variables 
and Hamiltonians used, we denote the independent variable as r. 
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Let h = h(z ) be an arbitrary function and z{r) be the solution of the 
system ([11]) with initial condition z(t ) = zq . It is clear that the derivative 
of the function h(r) = h(z(r)) for the Hamiltonian system ( |TTD has the 
form 

From ( |TB| ) we can obtain the expansion of h(r) with respect to a small 
At = t — t 

h(r) = h(r ) + Ar{h, H}(r ) + 0(Ar 2 ) (14) 



Writing /(r) and (7(r) instead /i(t) in ( p!4]) and calculating their Poisson 
brackets with respect to the variables z , we see that 

{/, g}(r) = {f, g}(r ) + 

+ At ({f, {g, H}}(r ) + {{/, H}, g}(r )) + 0(Ar 2 ) (15) 
Applying Jacobi identity to the terms of order At in flTBl) one finds 

{/, 9}(r) = {/, 9}(r ) + Ar{{/, g}, H}(r ) + 0(Ar 2 ) 
or, equivalently 

{/, 9}(r) - {/, sHto) 



At 



{{/, g}, H}(t ) + O(At) (16) 



The formula (|T6|) in the limit when At — > gives us the fundamental 
property of the solution of the Hamiltonian system: 

The flow of a Hamiltonian system preserves the Poisson bracket :Q 

-f- {/, 9} = {{/, 9}, H} (17) 



2 We use the term "preserve" to mean that the algebraic form is not changed. 
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2.3 Hamiltonian Extension of the Equations of Classi- 
cal Spin-Orbit Motion 

We now introduce the Poisson bracket^ 

{/(*)> 9(z)} = fq ■ 9p - fp ■ 9q + [fs x 9s] • s (18) 

in the 9-dimensional phase space z — (x, s) of 6 orbital variables x = (q,p) 
and 3 spin variables s and consider a Hamiltonian system of ordinary 
differential equations 

f - « ^ < 19 > 

with the Hamiltonian function 

H = H orbt (t,x) + H>(t,f) -s (20) 
In the variables q, p and s the system (|19|) can be written as 



dq = dHorU + d(W-s) 
dt dv dv 



dH orbt d(w-s) 



dt dq dq 

ds 



dt 



W xs 



(22) 



(23) 



and we will understand the equations (|21|)-(p3[) as the Hamiltonian exten- 
sion of the equations of classical spin-orbit motion ([!])- (@). 

Note that the matrix J(z) for the spin-orbit Poisson bracket flT8|) has 



3 Note that when applied to the spin variables only, the Poisson bracket ( |l8| ) gives the 
usual result { s,, Sj} — eijk Sfc. 
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the form 



J{z) 



( 











1 














o \ 
















1 































1 













-l 































-1 































-1 













































-s 2 






















S3 





Si 


V 




















s 2 


-Si 


o / 



(24) 



and is not a constant matrix (in contrast to the case of classical Poisson 
brackets), but depends on spin variables. The structural matrix ( |2~3| ) can 
also be written in the more compact block diagonal form 

J(z) = diag(J, J s {s)) 



where a 6 x 6 constant matrix J is the symplectic unit |12fl and 



JJs 





-s 3 

S'2 



S3 



-si 



2.4 Connection between the Triangular System and its 
Hamiltonian Extension 

We now wish to point out some ways for establishing the connections 
between properties and solutions of the system (PT])-(P^) and the initial tri- 
angular system ([]])- (|2|) (truncation procedures). We will do this without 
ascribing any physical sense to the spin dependent members in the 
right sides of equations (|2~1|), (|22]) []. 

Because the statements, which will be listed below, are connected not 
with the specific form of the spin-orbit Hamiltonian, but only with proper- 
ties of the spin-orbit Poisson bracket, let us consider an arbitrary smooth 
Hamiltonian function possibly depending nonlinearly on the variables s 



H 



H(t, z) 



H(t, 



(25) 



4 Note that some authors (see, for example |53| g4|) ascribe the spin dependent members 
on the right sides of equations (|2l])-(|2^) to a quasi-classical effect of the spin on the orbit 
motion, so that in that case no truncation procedures are needed. 
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and introduce, in correspondence to the Hamiltonian system 

JJs) erad ? if 



dx 
dr 



J gmd g H, 



ds 
dr 



the triangular truncated system defined as 



^ = J ■ (grad^fT)!^, 



ds 

= Js(s) ■ (grad s -#)| s -. =g 



(26) 



(27) 



a) If 



z(r, r , z ) = (x(t, r , x , s ), s(t, t , x , s )) 



is the solution of (|26| ) which passes through the point zq = (x , s*o) when 
r = t , then 



f*( r , r , z ) = (x*(t, r , f ), s*(r, r , f , s )) 



where 

£*(t, r , f ) 



£(t, r , z , 0) and s*(t, t , x , s ) = t^t 

OSq 



■ s 



s =0 



gives us the solution of 

b) If the system ( p6|) admits an invariant function V(r, z) which can 
be represented in the form 

V(r, z) = V m {r, z) + V> m (r, z) 

where is a homogeneous polynomial of degree m in variables s , and 



lim 

|sl-»0 |S' 



>m 







then V^j(t, z) is a first integral of the system fl2~7l) . 

c) If x(t, r , x ) == 0(t, t , x ) is a solution of the first of the equations 
(|27|) , then the system ( p7|) can be written as a family of Hamiltonian systems 
of the type ( p6|) depending on parameters (t , x ) with the Hamiltonian 
function 



W (t, 0(t, t , x ) \ ■ s where W(r, x) = (grad^i/) 



s=0 
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3 Hamiltonian Methods for the Extended Sys- 
tem 



3.1 Degenerate Poisson Brackets and Reduction of the 
Order of a Hamiltonian System 

If there are nontrivial functions (Casimir functions) fi(z) (maybe given 
locally on the manifold M) such that 



{fh h] 







(28) 



for any function h(z) then the matrix J(z) is degenerate and this Poisson 
bracket is said to be degenerate. (For a degenerate matrix J(z) of constant 
rank, the functions /; in (|2"B] ) locally always exist.) 

If all (at least all functionally independent) Casimir functions // have 
been found, then from properties ( |i3|) and (^8|) it follows that for an arbitrary 
Hamiltonian the trajectories of the system ( |TT| ) z(t) lie on the intersecting 
level surfaces 



fi(z 



Cl 



const (I = 1,2, ...,m) 



(29) 



where the Poisson bracket no longer remains degenerate. 

Locally any Poisson bracket of constant rank can be brought into the 
form (Darboux theorem) 



J 



( ° 



nn 
In 

ran 



In O ri 



(30) 



where Oki is a k x I zero matrix and I n is a n x n identity matrix, 
and dimM = 2n + m . Thus in local coordinates we obtain the classical 
Hamiltonian system with n degrees of freedom depending on m parameters 



l c l> c 2, 



Cm m 



]))■ 



3.2 Degenerate Poisson Brackets for Global Variables, 
or Local Darboux Coordinates ? 

The spin-orbit Poisson bracket flTBD is degenerate. It has the nontrivial 
Casimir function fi = \ s\ 2 and on the level surface f\ = const > its rank 
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is constant and is equal to 8. This means that we can decrease the dimensions 
of the system fl26|) by introducing Darboux coordinates. Thus we obtain the 
classical Hamiltonian system with 4 degrees of freedom depending 

— » 1 2 

on one parameter | s | . It is clear that the Darboux coordinates are not 
unique and may be introduced in various ways. We consider only one typical 
example. 

Let i, j, k be an arbitrary orthogonal system of unit vectors in three 
dimensional space R 3 satisfying the condition 

i ■ j x k =1 

We introduce three new spin variables ip, J, I by the equations 



s ■ i 



J 



< s-j = VI- J 2 cos(^) (31) 
s- k = y/I — J 2 sin('0) 

or, equivalently 

s[t/j, J, I) = J ■% + VI -J 2 (cos(^) • J + sin(^) • fc) 

Here J is the projection of the spin vector on the z-axis, ip is the polar 

I — * I 2 

angle in the transverse plane and I = | s \ .In the new variables the spin 
part of motion equations (|26|) becomes 

i> = Hj, j = -Hi,, I = (32) 
where the Hamiltonian (p25| ) takes on the form 

H(t, x, if;, J, 7) = H (r, f, s(ip, J, I)) 
Unfortunately, when 

((grad^)-/) 2 + ((grad s H)-k) 2 + (33) 

this coordinate system cannot be extended onto the whole sphere / = const > 
0, since it has a singularity for I — J 2 = 0. This means we need to have a 
whole atlas of local coordinates systems (at least two local coordinate systems 
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defined by different vectors i%, ji, k% and i 2 , j 2 , k 2 in (|3T|) ) for a complete 
description of spin motion on the sphere in the electric and magnetic fields 
depending on time and position of the particle. We will have the same diffi- 
culties with any other Darboux coordinates because they are defined by the 
topological properties of the sphere. So the way pointed by the Darboux 
theorem does not look like the most natural or straightforward ap- 
proach to the problem of investigation of polarized beam dynamics 
and we prefer to study the equations of spin-orbit motion using initial global 
variables x, s and the Poisson bracket (degenerate). 

Remark: Note that in Darboux coordinates considered here the Hamil- 
tonian (p0[) linear with respect to spin variables takes on the form 



H = H orbt + (W ■ i) J + ^T^P ((W-j) cos (VO + (W ■ k) sin ty) 
and the condition ([33]) now reads as 

(w-j) 2 + (w ■ k) 2 ^ 

3.3 The Properties of Solutions of the Extended Equa- 
tions of Spin-Orbit Motion Independent of the Spe- 
cific Choice of the Hamiltonian 

Which properties of solutions are independent of the specific choice of 
the Hamiltonian function in the equations of motion (|TT|) ? These are, for 
example, properties connected with the existence of Casimir functions and 
with the preservation of the Poisson bracket along the trajectories of the 
Hamiltonian system. 

The spin-orbit Poisson bracket ([18]) has the Casimir function | s | 2 . This 
means that the length of the vector s(t) is preserved during the motion, i.e. 

I s(r) | = | s(t ) I 



Let z(t) = 0(t, t , zq) be the solution of the system (|26|) , where 
0(tq, r , Zq) = z . Using the Taylor series expansion of the function 0(r, r , z ) 
with respect to spin variables we obtain: 

' x(t) = F(t, t , x ) + 0(| s |) 

(34) 

s(r) = A(t, r , x ) ■ s + 0{\ s | 2 ) 
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where A(t, t , xq) is a 3 x 3 matrix. The map (^) preserves the Poisson 
bracket (0). Using this property we find: 

a) The Jacobian matrix of the vector-function F(t, t , So) is symplectic: 



Here the 6x6 matrix J is the symplectic unit [|12|] and the symbol 'T' 
indicates transpose of a matrix. 

b) Every element of the matrix A(t, to, Sq) is equal to its own cofac- 
tor. For a 3 x 3 real nonsingular matrix this means that A(t, Tq, So) is an 
orthogonal matrix and det A(r, r , So) = 1 , i.e. 

A(t, t , So) g SO(3) 

We have briefly discussed the properties connected with the existence 
of Casimir functions and with the preservation of the Poisson bracket, but 
without doubt, Liouville's theorem on the conservation of volume is one of 
the most popular properties of classical Hamilton systems (at least in the 
accelerator physics) which is independent of the choice of a specific Hamilto- 
nian. Strictly speaking, this is a statement about the existence of an integral 
invariant (dim ^-dimensional) of the density g{z) = 1. 

It is well known JT3J that the system of ordinary differential equations 

f = f(T,z) 

with continuously differentiate right hand side admits a non-negative con- 
tinuously differentiable function g(z) (in particular the function g(z) = 1) 
which serves as the density appearing in an integral invariant if and only if 

div* (g(z)f(r, z)) = (35) 

For the Hamiltonian system ( fL^) the equality (|35D has the form 

&v z {g{z) J(z) gra,d g H(T,z)) = (36) 



It follows from (|36|) that the Hamilton system (12) will preserve the phase 



space volume for a given Hamiltonian function H if and only if 

^, m ^,H) = a A.^ = ^Ay^ = o (37) 
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Example: Defining in the two-dimensional Euclidean space the Poisson 
bracket by the equality 



{zi, z 2 } = z x 
we will have the Hamiltonian system 

OH dz 2 



dzx 
dr 



zx 



dz 2 



dr 



-Zx 



OH 

dzx 



Suppose H = a ■ z 2 . Then a solution of the system ([ 
of the formulae 

z 2 (r) = z 2 (0), zx(r) = Zl (0)-e° 
Consequently, an image of the single square 



(38) 

is given by means 



< 21 (0) < 1, 



< z 2 (0) < 1 



will have the area e aT after a displacement along trajectories. The area is 
preserved for a = , and — > for a < , and — > oo for a > 0. 

Examining the condition fl37|) we find that the Hamiltonian system will 
conserve phase space volume independently from the choice of Hamiltonian 
function if and only if 



(9 7 
^~dz~ 



0. 



J 



, dim z 



(39) 



For the spin-orbit Poisson bracket all the values ( |39D vanish and consequently, 
the phase volume in the 9-dimensional space is preserved. 

Moreover the phase volume is the integral invariant and in the reduced 
8-dimensional phase space, this is a direct product of the two-dimensional 
sphere of a fixed radius p > and the 6-dimensional space of orbital variables 
(level surface of the Casimir function | s*| 2 ). To proof this fact consider an 
arbitrary sufficiently small domain D and introduce the spherical coordinates 



s ■ i 



s " J 



p cos sin <p 
p cos 9 cos <p 



(40) 



s- k 



psin 9 
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where -§ < Q < §, 
satisfying the condition 



< if < 2n and where the unit vectors i, j and k , 



j x k 



1 and forming an orthogonal basis, are 
chosen in such a way that the coordinate transformation P0| ) is nonsingular 
within some open set which includes the domain D. In new variables the 
equations of motion (E6f) have the form 



^ = J grad^tf, 



dp 
dr 



-h = 



(41) 



dd 1 <9F dy> 

dr p cos 9 dip' dr 

If we neglect the equation 

d ~T = 

it is easy to check that the remaining system of eight equations has the 
integral invariant with the density 

g = p 2 cos9 (43) 

but the density (|^) is just the density of a volume in the 8-dimensional 
phase space considered. 

Remark: One can check that the properties listed above are correct 
not only for the extended but also for the original triangular system (in 
fact, this was the criterion for their selection). For the triangular system 
these properties can be regarded as properties independent of the choice of 
external electromagnetic field. 



dH 



p cos 9 89 



(42) 



4 Canonical Transformations and the Intro- 
duction of Machine Coordinates for Circu- 
lar Accelerators 

In the theory of circular accelerators it is useful to describe the spin-orbit 
motion in terms of a curvilinear coordinate system associated with the design 
orbit. In the previous sections we have introduced Hamiltonian extension of 
the equations of classical spin-orbit motion. Hamiltonian systems have a very 
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special form, and the special form is not preserved by an arbitrary change 
of variables. In this section we describe the transformations of phase space 
that are canonical with respect to Poisson bracket (|IBD and preserve that 
special form, and which allow us to make coordinate transformations using 
the Hamiltonian function ( p0| ) directly instead of the equations of motion. 

If we linearize the resulting Hamiltonian equations with respect to spin 
variables and then neglect the effect of spin on the orbit motion (triangular 
truncation procedure) we obtain the transformed version of the triangular 
system too. 



4.1 Canonical Transformations of Phase Space 



As was already noted, in a fixed local coordinate system z the Poisson 
bracket is completely defined if we know the values of the elements of the 
skew-symmetric matrix J = {{z i: Zj}) as functions of z. Before introducing 
the canonical transformations we will discuss the converse problem. Under 
what conditions on the skew-symmetric k x k matrix J(z) will the binary 
operation 

F j(f, 9) = grad^/ • J gmd^g 

be a Poisson bracket? This operation is automatically bilinear and antisym- 
metric and satisfies the Leibnitz rule. Hence the operation Fj(*, *) will be 
the Poisson bracket if and only if it satisfies the Jacobi identity. The Jacobi 
identity written in terms of elements of the matrix J has the form (see, for 
example |T5| |) 



J, 



ml 



dJij j dJ m i ~ dJj m 
+ Jji ■ — — + Ju ■ 



■„ ., -„ ., 0, i,j,m = 1, ...,k (44) 

dyi dyi dyi 

Returning to the main theme, consider a coordinate transformation from 
old variables z to new variables y in the Hamiltonian system of ordinary 
differential equations fll2D 

z = 0(y) (45) 
with a nondegenerate Jacobian matrix (dim z = dim y = dim M = k) 

I dipi _ _ _ dipi \ 
dyi dy k 



\ dyi 



dyk ' 
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If we take into account the connection between the gradients of the Hamil- 
tonian function H in old and new variables 



' d(p\ T 



grad^ H(t, (p{y)) = grad z -if(r, z 



(46) 



z=0(y) 

we easily obtain from ([12]) the differential equations for y 



dr 



We now introduce the function H(t, y) = H(t, 0{y)) and the skew- 
symmetric matrix 

and rewrite the system (fl"7D in the form 

^ = J(y)gvad g H (48) 

Equation (|48|) resembles a Hamiltonian system with a new Hamiltonian func- 
tion H(r,y). But in fact it will actually only be a Hamiltonian system (in- 
dependently from the choice of specific Hamiltonian) if the matrix J satisfies 
the Jacobi identity (f4"4"|). This condition looks very complicated, but there 
are two important situations when it becomes trivial: 

1. The matrix J(y) is constant (independent from y). 

2. The matrix J(y) is equal to the matrix J(y)- 

Example 1: For the Darboux coordinates ([H|) the matrix J has the 
form ( pUD and hence is constant. 

Example 2: Introduce the new spin-orbit variables by the equations 

•Cold F i^new) 

(49) 

Sold Snew 

If the Jacobian matrix of the vector-function F is symplectic then the matrix 
J is equal to the matrix J (E3). 
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Now we are ready to introduce the necessary definitions. 

Definition 1: A coordinate transformation z = (p{y) is a Poisson 
transformation of phase space if the matrix J{y) satishes the Jacobi 
identity (HJ. 

Definition 2: A coordinate transformation z = (p{y) is a Canonical 
transformation of phase space if the matrix J(y) is equal to the matrix 

•/(//)• 

The Poisson transformations preserve the Hamiltonian form of the initial 
system, but their subset, canonical transformations, has additional helpful 
properties. Rewrite the condition J(y) = J(y) in the form 

*«» ■ (I) j ^ (IT <»> 

The equality (|50|) means that the image of the Poisson brackets of coor- 
dinates functions is equal to the Poisson brackets of the images. Thus a 
canonical transformation is an "isometry" of phase space, and we can find 
the Hamiltonian equations in new variables by just making the coordinate 
transformation in the Hamiltonian function (as in the usual way). 

Example 3: The solution z(t) = <p(r, To, Zo) of the Hamiltonian system 



11) for any fixed r > tq is a canonical transformation of phase space: 



A<P(r,r ,z )) = J(z ) (51) 



Remark 1: The condition for the map ( f45D to be symplectic is often 
used not in the form following from (]5"0|) 




= J (52) 

where the matrix J is the symplectic unit, but in the form 

( 9<p s 




= J (53) 
It can be easily shown that since J 2 = —I these two conditions are equiv- 



alent, but it is not true in general that ( j52|) and ([53]) are equivalent even for 
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a constant nondegenerate matrix J in (|5T)|). To illustrate this let us intro- 
duce a Poisson bracket in four- dimensional Euclidean space by means of the 
structural matrix 



J 



and consider the map 
Zi = 2/1, z 2 



( 





i 









-i 






















a 


\ 








—a 


/ 



const 



2/2 - 2/4, 



^3 



a 2/i + 2/3, 24 



2/4 



which is generated as a displacement along trajectories of the Hamiltonian 
system with H — z\ • Z4 for the time r = 1 . Calculating the Jacobian matrix 

) gives 



of this map we can check that fl52|) is satisfied but that 

/ 




J + 



a 2 -1\ 
OOOO' 

V 1 - a 2 / 



^ J if 



^ 1 



Remark 2: Let a Poisson transformation satisfy the equality 

JiM)) c[^ \ .nil 



(54) 



for some nonzero constant c. If the relation ( p4[) holds, then in the new local 
coordinates the structural matrix of the Poisson bracket is determined by 
the matrix J(y) = c J(y) ■ If c ^ 1 this transformation is not canonical. 
Sometimes for c / 1 it is useful to introduce for the transformed system a 
new Poisson bracket, determined by the matrix J(y) = J{y) ■ Defining 
the new Hamiltonian by means of the rule 

H(t, y) = c-H(t, (p{y)) 



we can consider the transformed system as being Hamiltonian with the 
"same" Poisson bracket. For instance, in accelerator physics this method 
is used in studying betatron oscillations when the transverse momenta are 
normalized by the value of the kinetic momentum of a reference particle, 
which is assumed be a constant. 
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4.2 Nonautonomous Canonical Transformations 



Except for transformations of the type (|45|), in this paper more general 
transformations of variables are used 

z = <p(r, y) (55) 

depending on r as a parameter. Then we need to answer the question: "In 
what case is the transformation (|55|) canonical and how is the new Hamilto- 
nian to be calculated?" 

Definition 3: A coordinate transformation z = (p(r,y) is a nonau- 
tonomous canonical transformation if for any fixed r the condition 
(|5^ ) holds and there is a differentiable function F(r, y) such that 

^.^^-g;)" 1 .^ (5 6 ) 

Then the new Hamiltonian is given by the formula 

H{t, y) = H(t, <p(r, y)) + F(t, y) 



Example 4: Defining in the two-dimensional Euclidean space the Poisson 
bracket by the equality 

{zi, z 2 } = z 1 + z 2 
it easy to check that the coordinate transformation 

Z\ = (1 + b)yi + by 2 + c 

(57) 

, Z2 = -byi + (1 - b) y 2 - c 

is canonical independently from the choice of constants b and c. Letting now 
b and c be functions of r we find that for the set yi + y 2 ^ the equation 
(|56| ) has the solution 

F(r,y 1 ,y 2 ) = -(yi + y 2 )-^ - sign(yi + y 2 ) In \y x + y 2 \ ^ 

and that this solution is unique up to an additive arbitrary function of the 
variable r. This means that if 
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we cannot consider the transformation fl5"T| ) to be a nonautonomous canonical 
transformation because the equation (|56|) does not have different iable solu- 
tions (at least in the neighbourhood of the set yi + yi — 0). So for general 
Poisson brackets the solvability of the equation (|56|) does not follow from 
satisfying the condition (|50"D for all values of r. 

We will not set ourselves the target of studying the general properties 
of nonautonomous canonical transformations, but instead we consider some 
examples directly connected with the purpose of this paper. 

Example 5: Let z[r) = 0(r, tq,zq) be a solution of a canonical system 
with the Hamiltonian F(t, z) , so that 



q- = J(<f>) -grad £ F(r, z 

Taking into account ( 46]) and ( |5lD one obtains from ( [58] ) 
-l 



(58) 



_ I = J(%o) "grad^ F(t, ${t, r , z )) (59) 



From (|59|) and (^l|) it follows that a coordinate transformation 

z = 0(r, r , y) 

is nonautonomous canonical and the new Hamiltonian is given by 

H(t, t , y) = H(t, 0(r, r , y )) - F(r, 0(r, r , y )) (60) 

As a particular case, the formula (^) contains the Hamiltonian version of 
the method of variation of constants when H = Hi + H 2 and the function 
F is chosen to be Hi . 

Example 6: Generalizing example 2 consider the transformation: 

•^old •Knew) 
Sold $new 

With the assumption that the Jacobian matrix of the vector function F(r, x) 
with respect to the variables x is symplectic for all values of r , the equation 
(|56|) is reduced to 




grad^F = — J — — (61) 
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The same symplecticity condition allows us to show that the matrix 

/ 



d_ 

dx 



dF 



OF 



\ 



dx 



dr 



is symmetric and, consequently, that the equation fl6l]) has a solution (defined 
up to some additive Casimir function). 

Omitting a proof, we also point out that if only orbital variables are 
transformed, then we can use the classical technique of generating functions. 

Example 7: The linear transformation of spin variables 




x t 



A{r) s r 



(62) 



will satisfy fl50|) if and only if A e SO (3) for all values of r. The condition 
(|56| ) now becomes: 



V 






Ml 


dF \ 

ds 2 


ds 3 


dF 





dF 


ds 3 


ds i 


dF 


dF 


o ) 


&S2 


dsi 



A 1 



dA 



(63) 



where the function F does not depend on x. 

Taking the derivative with respect to r in the identity 



A\r) A{t) 



we find that the matrix 



A 



dA 
dr 



is skewsymmetric, and hence that (|63| ) has a solution which can be expressed 
as follows 



dr J 



-dA s 
dr 



A T ^ 



13 



■ S 3 



12 



So ( |62"D will be a nonautonomous canonical transformation for an arbitrary 
differentiable matrix A G SO (3). 
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4.3 The Coordinate Frame Connected with the Closed 
Design Orbit 

In this section the words "the closed design orbit" mean some suitable 
closed curve which has a continuous unit tangent vector. Let the closed 
design orbit be described by the vector fo(z) , where z is the length along 
this curve. Supplement the unit tangent vector 

dz 



T 



with two unit vectors N and B satisfying the conditions 



B = 
or, equivalently 



T x N 



N = 



TxN 



BxT 



■ B 



T = 



N x B 



The triplet T , N , B thus forms an orthogonal right handed coordinate 
system. We assume that the evolution of this coordinate system as the 
variable z changes is described by a periodic solution of the system of ordinary 
differential equations of the Fresnet type 



dT 
dz 
dN 
dz 
dB_ 
dz 



— = -hN-aB 



— = +hT + 3iB 



+aT - aeiV 



We will find the transformation from the old spin-orbit coordinates to the new 
spin-orbit coordinates connected with the vectors T, N , B as a composition 
of two successive transformations: the first one changes the orbital variables 
and second one changes the spin variables and the longitudinal momentum. 

4.3.1 Transformation of Orbital Variables 

In the new coordinate system an arbitrary orbit-vector f lying in a 
sufficiently small neighbourhood of the closed design orbit can be written in 
the form 



f (z) + xN + yB 
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We will take the parameters z, x, y to be the new orbital variables. The 
transition from the old coordinates to the new coordinates is made, as usual, 
with the help of the generating function depending on the new position and 
old momentum variables 



F(r, p ) = — \To(z) + x N + y B J ■ p 

The new momenta are given by the equations 
OF 



dz 



p- (Cl + hx + ay)f + ae (xB - yN)) (64) 



OF 

P x = = P- N (65) 

dF 

P y = -— = p-B (66) 

This transformation is a canonical transformation with respect to the spin- 
orbit Poisson bracket. 



4.3.2 Transformation of Spin Variables and Longitudinal Momen- 
tum 

The new spin variables are introduced via the equation 

Sold = C(z) s new (67) 

where C(z) = \N(z), B(z), T(z)^j is a 3x3 matrix. 

Since the matrix C depends on the variable z the coordinate transforma- 
tion (67]) is not canonical. To make it canonical we change the longitudinal 
momentum too: 

pold = pnew _ + hSy _ (gg) 

Here s x , s y , s z are the components of the spin vector s new . 

The coordinate transformation fl67D , (|6S|) satisfies the condition (0) and 
hence is canonical. Since the old Hamiltonian contains all vectors as scalar 
and cross products then the new Hamiltonian will have the same form as 
the old one if we imagine that all vectors in (|20| ) are written in terms of 
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the projections on the unit vectors N, B, T and we take into account the 
formulae 



Pn = P - N = P x 



Pb = P- B = P v 



Pf = p-T 



1 



1 + hx + ay 



(P z - ee(xP y - yP x ) - a s x + h s y - £es z ) 



which we can easily obtain from (Bl)-(pBD and (55). Here we have reverted 
to using the symbol P z instead of p™ ew . 



4.4 Change of Independent Variable in Nonautonomous 
Hamiltonian Equations 

Let the right hand part of the first equation of the Hamiltonian system 



TT| ) satisfy the condition 



dr 



{z u H(t, z)} ± 



(69) 



This means that the variable Z\ (r) changes monotonically with changing r 
(strictly increasing or strictly decreasing) so that one can introduce it as new 
independent variable. In many cases of practical importance the new "time" 
scale connected with Z\ gives us certain advantages and we wish to discuss 
the procedure of its introduction for Hamiltonian systems. 
For convenience we introduce following notation 



Z = f (q, P, Vi, V2, - Vk-2) = (q, P, y) 



(70) 



for the components of the vector z and assume that the matrix J for the 
Hamiltonian equations (0) has the form 



J(z] 



( 1 
-10 


^00 



J(v) 



\ 




(71) 



This is not a restriction, because on the one hand the spin-orbit Poisson 
bracket (|18[) has the necessary form and on the other hand any matrix J can 
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be brought into this form if the condition 
the system (|TT|) becomes 



holds. With the new notations 



dq dH dp dH dy, t - /^\^H 
dr dp dr dq' dr 13 dyj ' 



k-2 



(72) 



and the condition fl69|) now reads as 

dH 

dp 



^ 



(73) 



We will interpret the procedure of changing the independent variable as a 
procedure of the reduction of an autonomous Hamiltonian system to a family 
of nonautonomous Hamiltonian equations of smaller dimension defined on the 
level surfaces of the initial Hamiltonian function. With this aim in mind we 
introduce two additional canonical variables (E, v) and the new Hamiltonian 



H{v, E, z) = H{v, z\ 



E 



to obtain an autonomous Hamiltonian system in a higher dimensional phase 
space 

dq dH dp dH 

dr dp ' dr dq 



(74) 



dyj 
dr 



My, 



dH 



k-2 



(75) 



dE dH 

dr dv ' 



dv 
dr 



dH 
"dE 



(76) 



Using the condition (|73|) we obtain the differential equations for the new 
independent variable q 



dyi, dy^ dr - dH/dyj 

dq dr ' dq tAV) dH/dp' % 



l,...,k-2 



(77) 



dE dE dr dH/dv 
dq dr dq dH/dp 



dv dv dr dH/dE 



dq dr dq 



dH/dp 



(78) 
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We now wish to show that the equations ([77]), ( |7S|) are the family of nonau- 
tonomous Hamiltonian systems defined on the level surfaces 



H{E, v, q, p, y ] 



Co 



(79) 



Let 



P = K(E, v, y, q, c ) 

be the solution of the equation ([fj|). According to the implicit function 
theorem this solution exists if the condition ([FJ|) holds. Taking the derivative 
with respect to the variable yi in the identity 



we have 



H(E, v, q, K(E, v, y, q, c ), y) = c 



m on dK 

dyi dp dyi 



(80) 



(81) 



From (81) it follows that 



Similarly we obtain 
dK 



dE 



dK 

dyi 



8H/dE 

m/dp '■ 



dH/dy, 
dH/dp 



dK 

dv 



(82) 



dH/dv 
dH/dp 



(83) 



Comparing (|77|) , (|78|) and (|82D , (|83D and introducing the Hamiltonian func- 
tion 



H 



K{E, v, y, q, c ) 



we see that the equations (|77|) , (|78|) become the family of nonautonomous 
Hamiltonian systems depending on the parameter cq 



dE _ dH dv 
dq dv ' dq 



dH dyi - ,-,xdH 
dE' lq~ = Jij( - y) dy~' % 



k-2 



J4) 
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After solving the system (Q) for a fixed value of Cq we can find the dependence 
of p on q using the identity 



p(q) = K(E(q), v(q), y(q), q, c ) 
and then determine q as a function of r from the equation 

r~r = J^L (85 ) 

where 

&J~l 

(equation fl35|) follows from the first of equations (ff^)). 

Now we have to remember that we did not start from an autonomous 
system (fF4])-(|76|) but from nonautonomous system (|72|). This means that 



we do not need equation fl85[) because we have the dependence r(q) from 
the second equation of fl84|) {v is just another notation for r). We also 
have freedom in the choice of the initial condition for the variable E. This 
means that we can choose a single fixed value of the parameter Co (usually 
Co = 0) and replace the initial nonautonomous Hamiltonian system by the 
nonautonomous Hamiltonian system flS^ ) with the same matrix J as in ([71]) . 

4.5 Length Along the Design Orbit as Independent 
Variable 

For a circular accelerator the spin-orbit Hamiltonian is always a periodic 
function of the variable z, but its dependence on time may be more compli- 
cated (for example, in the acceleration mode). This is one of the reasons for 
introducing the length along the design orbit as the independent variable. 
We will assume that 

dH , 

m * < 86) 

The condition (^) approximately means that during the motion in the ac- 
celerator the particle cannot reverse its direction. 
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Following the previous subsection we introduce two additional canonical 
variables (E, r)Q and the new Hamiltonian 

T~L H(l~, X, P x , y, Py, Zj Pzi $xi Syi Sz) E 



dE dH dr &H _ 

~dt ~ ~dr' ~dl ~ ~~dE 

Now we need to solve the equation 

H{E, r, x, P x , y, P y , z, P z , s x , s y , s z ) = (87) 
with respect to variable P z to obtain the new Hamiltonian 

H = -P Z (E, t, x, P x , y, P y , s x , s y , s z , z) (88) 

The dependence of the Hamiltonian 7i on the variable P z is more complicated 
than in the pure orbital case, but nevertheless, we can solve the equation (p7| ) 
with any required precision with respect to spin variables using the method 
of successive iterations. 



4.6 The Hamiltonian in New Variables up to First Or- 
der with Respect to Spin Variables 

In this subsection we discuss the general form of the spin-orbit Hamil- 
tonian in the new variables (E, t, x, P x , y, P y )Q up to the first order with 
respect to spin variables. This Hamiltonian is 

H = H orbt (E, t, x, P x , y, P y , z) + H spin (E, t, x, P x , y, P y , z, s) (89) 

with 

H orbt = -ee(xP y - yP x ) - (1 + hx + ay)- 



5 Reserving for a while the new symbol r for the time t appearing in the spin-orbit 
Hamiltonian, and still reserving the symbol t for the independent variable. 

6 Here and in the Hamiltonian (^8|) we have reverted to using the symbol t instead of r. 
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8ex7Tg + asy-Rtf — (1 + hx + ay)nf — - A z 



and 



H, 



spin 



— a s x + hsy — as s z + 



+ 



[1 + hx + ay){E - e$) 



W-s 



C \/ r 2 ~ TUqC 



- ~ c A^j - (Py - ^Ag 



m 7 



as x + hs y — ?BS Z + (1 + hx + ay) • s 



7Tj; 



Here W has the same form as in (20) with if, B, S written in terms of 
projections on the vectors N, B, T 



71 



N 



P — - A - 

1 x /V> 
C 



71 



B 



7X n 



(E-e$) 



rule- - [P x --Aff) - [Py - -A a 



'(E-e$Y 



1/2 



2 2 2 2 

m C ~ ~ n B 



and we use the notation 

A z = (1 + hx + ay) Af + se(xAs — y A^) 
The value of 7 is defined now through the new canonical variable E 



7 



E - e$ 
m c 2 



To obtain the Hamiltonian ( p9|) we have used the condition (86) in the form 

dH 



dP z 



> 
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so that the vector T is chosen with the same orientation as the direction of 
particle flight in the accelerator. 

To complete the description we also give the expressions for the projec- 
tions on N, B, T of the electric and magnetic fields in terms of the vector 
and scalar potentials. 

The magnetic field: 

N ~ 1 + hx + ay { dy dz * y ° T ) 

B s = —± (?£L-?£- +tBxBi ) 

1 + hx + ay \ dz ox J 

The electric field: 



£s = 



dx 


dy 


d$ 


1 dAfi 


dx 


c dt 


.9$ 


ldA s 


dy 


c dt 




( 9$ 


+ ae 


\ dx 



1 (m ( 9$ d$\\ 1 dA f 
^ 1 + hx + ay \dz 86 \ dx X dy J J c dt 

(Useful formulae which allow us to get the equations for vector and scalar 
potentials in our curvilinear coordinate system may be found in the Appendix 
A). 



5 Linear Differential Equations of Spin Mo- 
tion 

In this section we discuss the situation in which the Hamiltonian function 
is linear in the spin variables and does not depend on the orbit variables. In 
particular, this case includes the description of the behaviour of the spin 
vector on the closed (or any other chosen) trajectory of orbital motion, using 
the triangular system. 
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5.1 Matrix Representation of the Hamiltonian Func- 
tion and Simple Properties of Solutions 

The Hamiltonian function, which is linear in spin variables and does not 
depend on orbit variables, has the form 

H(t, s) = w(t) ■ s (90) 

We now construct the skewsymmetric matrix C(w) according to the rule 

/ -w 3 w 2 \ 
C(w) = \ w 3 — w\ 

\ -W 2 W-l / 

where w±, w 2 , W3 are the components of the vector w 

w = (wi, w 2 , w 3 ) 

It is easy to verify that with the help of the matrix C(w) the Hamiltonian 
(p0| ) can be written as follows 

H(t, s) = is -curl s (C(w) ■ s) (91) 

and the equations of motion take the form 

^ = C(w(t)).s (92) 

The matrix representation ([Jl|) of the Hamiltonian function (0) resem- 
bles the representation of the Hamiltonian of linear orbit motion as a quadratic 
form and is particularly convenient for linear nonautonomous canonical trans- 
formations of variables. Introduce, for example, the new spin variables u by 
the equation 

s = A{t) u, A{t) e SO(3) (93) 

Substituting (^) in the equation of motion (|9~2"D, we obtain 

du 
d^ 



AC A - A 



■ dA s 
dr 



u 
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where 

A T C(w)A - A T ^ 
dr 

is a skewsymmetric matrix again. Consequently, the Hamiltonian function 
in the new variables u can be written in the form 

H(r,u) = ~ u ■ curl £ ^ ^4 T C A - A T< ^j ■ uj (94) 
The following simple properties of the matrix notation are almost obvious: 



a) curl s - ( A ■ s ) + curl s - ( B ■ s ) = curl s - ( (A + B) ■ s) 

b) A T C{w) A = C(A T -w) for A e SO(3) 

Let M(r, ro) be the fundamental matrix solution of (p2]). It has been 
mentioned in subsection 3.3 that for any r > r 

M(r, r ) G SO(3) (95) 

If the vector w is constant (does not depend on r) and \w\ ^ , then the 
matrix M is determined via the matrix C(w) from the formula 

M(t, r ) = exp((r-r )C(«;)) = 



J + sin(| w\(t- r )) • + (1 - cos(| w \ (r - r ))) • f (96) 

| it; | y | w | y 

In the case \w\ — the matrix M is equal to the identity matrix. 

The properties of the solutions of the system of linear differential equa- 



tions (|92|) given below follow from the Hamiltonian character of the system 
(l92f) (in the case studied here it is equivalent to the fact that the matrix 
C(w) in ( p2|) is real skewsymmetric): 

a) The norm of every solution is conserved. 

b) The angle between any two solutions is conserved. 

c) If n{r) and m(r) are solutions of the system (0), then 
/(t) = [n{r) x m(r)] is a solution too. 
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d) The function n(r) ■ s is a constant of motion of the system (£55J) 
if and only if n(r) is a solution. 

e) If the vector w(t) in ( |90D has period T in r, then the system fl92|) 
has a T-periodic solution. 

f) If the vector w{t) in fl90[ ) has period T in r and the system (|92| ) 
has two linearly independent T-periodic solutions, then any solution of the 
system (|92"D is T-periodic. 

5.2 Linear Equations of Spin Motion with Periodic Co- 
efficients 

Now we consider a Hamiltonian function (0) in which the vector w has 
period 2ir in r, that is 

w(t + 2n) = w{t) 

We shall attempt to obtain a Hamiltonian function (|94|) of the simplest form 
by means of the linear coordinate substitution ( |93| ) with a 27r-periodic matrix 
A{t). 

5.2.1 Some Properties of SO(3) Matrices 

Recall that the symbol SO (3) denotes the group of 3 x 3 real orthogonal 
matrices with the determinant equal to 1. 

Let A G SO(3). Since A' 1 = A T , the spectrum of this matrix is symmet- 
ric about the unit circle (that is, if \i is an eigenvalue of a SO (3) matrix, then 
so is /i -1 ). From the fact that the matrix is real it follows that the spectrum 
is also symmetric about the real axes. Taking into account the condition 
detv4 = 1, we obtain that the eigenvalues of the matrix A have the form 

\i\ = 1, [12,3 = cos A ± i sin A 

where 

cos A = - (Tr(v4) - 1) 

The next three lemmas give us some additional helpful properties of SO (3) 
matrices. 
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Lemma 1: A 3 x 3 nonsingular real matrix A belongs to the S0(3) group 
if and only if every element of the matrix A is equal to its own cofactor. 

Note, if every element of a 3 x 3 real singular matrix A is equal to its own 
cofactor then A = 0. 

Lemma 2: All the eigenvalues of a matrix A e SO(3) are distinct if and 
only if A ^ A T (or, equivalently, A 2 ^ I). 

Lemma 3: If all the eigenvalues of a matrix A e SO (3) are distinct then 
the nonzero vector 

— * 

n = -=r, k = (a 32 - a 23 , a 13 - a 31 , a 2 i - a 12 ) 
\k\ 

is a unit eigenvector of the matrix A corresponding to the unit eigenvalue, 
i.e. An = n. 



5.2.2 Real Jordan Canonical Forms of SO(3) Matrices 

Let n be an eigenvector of A corresponding to the unit eigenvalue 



An = n, 



n I = 1 



Supplement the vector n with two unit vectors rh and I satisfying the 
condition 



m 



I x n 



to form an orthogonal basis. Then the matrix 

B = (fh,l,n) e SO(3) 

— * 

In the new basis, constructed with vectors rh, I and n , the matrix A will 
have the form 

/ Am ■ rh Al • rh \ 
Am- 1 Al-l 
1 



A = B~ l AB = B 1 AB = 



\ 



e SO(3) 



It follows from the condition A e SO (3) that for some angle ip we can repre- 
sent this matrix in the form 

/ +cos(-^) +sin('0) \ 
A — — sin(^) + cos(tp) 

V o o i J 
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Using the equality AB = BA , we obtain 

A ■ m = cos(V') ■ m — sin(^) • I 

A ■ I = sin(-?/>) • m + cos(t/>) • / 

or, in complex form, 

A ■ m + i I ) = exp(i ip) ■ ( rh + il^j (97) 

It follows from (|97|) that exp(±z^) are the eigenvalues and m ± il 
are the corresponding eigenvectors of the matrix A (in particular, ip = 
±A (mod 2vr)) 

The matrix A is a real Jordan canonical form of the matrix A. We 
will denote it by A\. All possible canonical forms are given by the list 



f cos(^) +sin(?/>) 
,-T) = [ -sin(^) +cos(^) | . A-> 
1 

+ cos(^) + sm(ip) 
A-a = \ 1 I . .4 

— sin( , 0) + cos("0) 

1 
A- = ! +cos(^) +sin(V>) ! . A r , 
— sin^) + cos(tp) 



I + cos(^) — sin(^) 
+ sin(^) + cos(^) 

V o oi 

/ + cos('0) — sin(^) 

1 
\ + sin('0) + cos(V') 

/ 1 

+ cos(^) — sin( , 0) 
\ + sin(^) + cos(V') 



It is clear that all matrices A\ are similar one to another. For example, 



BjABi 

where Bi G SO (3) and 

10 
B s = | 1 | . B 2 
1 



A; 



6 



/ i 
o 

\ 



B, 





-1 






-1 





B* 



B* = 







-1 



We note that if ip — (mod27r), then all matrices A{ are equal one 
to another (and equal to the identity matrix). If ip — 7T (mod 2tc), then 
A l = A 2 , A 3 = A A and A 5 = A 6 . 
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5.2.3 Skew-symmetric Real Logarithm of SO(3) Matrices 

If the determinant of the matrix A is nonzero, then 

ln(A) = K 
is defined solution of the equation 

exp(K) = A (98) 

We now consider the problem of finding skewsymmetric real solutions of 
equation (|9~8|) in the case, when A G SO (3). 

Lemma 4: If A G SO (3) and A ^ A T , then all real logarithms of the 
matrix A are skewsymmetric matrices and can be expressed by formula 

HA) = ^ (a - 

where A and X* are complex conjugate eigenvalues of the matrix A distinct 
from unity. 

If A = A T then not all real solutions of equation (|9~8"D are skewsymmetric 
matrices. For example 

/ wd \ f -I 00\ 
exp -w/d 0= -1 e SO(3) 
V 00/ V 1 J 

for all real d ^ and all real w = 7r(mod27r). 

Lemma 5: Let n, m, and I be the orthogonal basis connected with 
matrix A G SO (3) and defined in the previous subsubsection. Then all real 
skewsymmetric logarithms of the matrix A are given by the formula 

I -n 3 n 2 \ 
\n{A) = {uj + 2nk) • n 3 -n x \ = (u + 2nk) ■ C{n) 

\ -n 2 rii / 

where real u (0 < uj < 2n) satisfies the equation [] 



exp\iu) 



- ( m + % f) ■ A ( m + i f) (99) 



7 Here and further on for complex vectors u, v £ C n we define u-v — Ui-v*+. . .+M n -w*. 
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and k is an arbitrary whole number. 

Note that the right hand side in ( |99D is just another expression for one of 
the eigenvalues (which we say is "conjugate with vector n") of the matrix A, 
and hence the equation ( |99|) always has a real solution. 

Below in this paper we will use the notation ln s A for real skewsymmetric 
values of the function In A. 

Remark: In accordance with Lemma 5 

ln s / = Ink ■ C{n) 

where n and k are an arbitrary unit vector and an arbitrary integer respec- 
tively 

5.2.4 Normal Forms for Hamiltonians of Linear Equations of Spin 
Motion with Periodic Coefficients 

For simplicity fix r = and let M(r) be the fundamental matrix of 
solution of the system (p2|). Then 

M(t + 2tt) = M(t) • M(2vr) 

Definition 1: A transformation of variables 

s = A(t)u, A{t) e S0(3) 

is called a normalizing transformation, if 

a) A(t + 2vr) = A(t) 

b) the Hamiltonian in the variables u does not depend on r 

c) the matrix A T (2ir) ■ M(2vr) • A(2tt) is a real Jordan canonical form 

Definition 2: The result of the application of the normalizing transfor- 
mation to the Hamiltonian with the periodic vector w(t) is called the 
normal form of the Hamiltonian for linear equations of spin motion with 
periodic coefficients. 

Let M be one of the real Jordan canonical forms of the matrix M(27r) 
and let the matrix B G SO (3) be such that 

B t M(2it)B = M 
41 



Lemma 6: The matrix 



A(t) = M(t) exp [ -— ln s M(27r) ] B 

V 27T / 

defines the normalizing transformation, which takes the initial Hamiltonian 
(9C) into the normal form 

H{$) = ^s-cmu(^(ln s M)-i) (100) 

Examining all possible real Jordan canonical forms of the matrix M(2tt) 
and explicitly calculating the logarithm in (|100|) we get 

Lemma 7: The Hamiltonian with the 2w -periodic vector w(t) can 
be reduced to the following normal forms 

H(s) = (±A + k) ■ s m (101) 

where 



A = — arccos -(Tr(M(2n)) - V 
2ir \2 

and the sign '+' or ' — in front of X, m 6 {1, 2, 3} and the integer k can be 
chosen arbitrarily. 

The number ±A + k in ( |101[ ) is called the spin tune. Using the freedom 
of the choice of k and of the sign ± in front of A we can normalize the spin 
tune to a value lying in the region from to 0.5. 



Note, that if A ^ I then the set of Hamiltonians ( |101| ) contains all possible 
normal forms. This is easy to see from the fact that if the Hamiltonian 

^ s ■ curl g (C ■ s) 

is a normal form, then the matrix exp(27rC) is a real Jordan canonical form 
of the matrix M(2ir). 

Another way to bring a 27r-periodic Hamiltonian linear in spin to the 
simplest form with help of a 27r-periodic coordinate transformation can be 
extracted from the results of Appendix B. 
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5.3 Connection between the SO (3) and SU(2) Groups 
and the Lax Form of the Equations of Spin Motion 

The fundamental matrix solution M(r, r ) of the linear differential equa- 
tions of spin motion fl92| ) is a 3 x 3 matrix and consists of 9 elements, but 
since M(r, r ) G S0(3) it can be described completely with the help of a 
smaller number of parameters using the connection between the SO (3) and 
SU(2) groups. Usually in accelerator physics this connection is established 
in the framework of the spinor formalism |16], [17|]. The method we shall use 
in this subsection (following ||) gives us the same results, but does not use 
the concept of spinors (at least in explicit form). 

Let us recall the definition of the SU(2) group: 

a 2 x 2 matrix U with complex coefficients belongs to the SU(2) group if 
detU = 1 and U ■ U* = I (the asterisk V indicates complex conjugation 
of a matrix). 

From this definition it follows that any matrix U G SU(2) has the form 

U = ( £ V a- a* + b-b* = 1 (102) 

We now define the matrix 

L = ( S \ Sl + iS * ), L* = L (103) 

\ S%- IS 2 -S 3 / 

corresponding to the vector s and introduce the ant i- Hermit ian matrix B 

B= % -.{ W \ ^ + B * = -B (104) 

2 \wt- iw 2 -w 3 J 1 

By means of the matrices L and B one can write the equations of spin motion 
(|92| ) in the form of a Lax equation 

^ = B ■ L - L- B (105) 

Note that the right hand side of equation ( 105|) satisfies the condition 
(B ■ L — L ■ B)* = B ■ L — L ■ B 
Lemma 8: If the matrix U(r, r ) satisfies the equation 

^ = B-U, U(r , r ) = I (or U(r , r ) = -/) (106) 
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then U(t, t ) G SU(2) and the solution of ( |i05| ) is given by the formula 

L(r, r ) = U(t, r ) ■ L(r Q , r ) ■ 17* (r, r ) (107) 



Remark 1: Because ( |1U7| ) is a similarity transformation it follows that 
for any r > tq the matrices L(r, tq) and L(ro, To) have the same eigenvalues. 
In our case this means that for any r > r 



ST 



s(r ) 



Remark 2: The equation of spin motion fl92|) can be written in the form 
(|105|) using another choice of the matrices L and B, different from ( |103|) and 
( |104j ). See, for example, 0. 

The following formulae give us the connection between the matrix U(t, tq) 
in ( |107| ) and the fundamental matrix solution of ( |92"D M(r, tq): 



.1/ 



/ Re(a 2 -b 2 ) -Im(a 2 + b 2 ) -2 ■ Re(ab) 

-2 • Im(ab) 



Im(a 2 — b 2 ) 
V 2 • i?e(a6*) 



i?e(a 2 + b 2 ) 
-2 ■ Im(ab*) 



aa 



bb* 



Introduce the quantities 
1 



m n + rn 22 + % (m 2 i - m 12 )) 



(108) 



1 



w 



(m 22 - mn - i{m 2 i + m 12 )) 



where m^i are the elements of the matrix M. Since for the orthogonal matrix 
M we have 



w ■ w + v ■ v 



^(m^ + m\ 2 + m\ x + m\ 2 ) ^ 



then using (|102|) and ( |108| ) we can define the elements of the matrix U by 
the equations 



v , b 



jr(m 13 + im 23 ), if v ^ 0, 



b = ±-y/w, a = — ^(mi 3 + z 777.23), otherwise. 
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Remark 3: The freedom of the choice of sign in these formulae (and also 
in ( p.06| )) is connected with the fact that the SU(2) group overlaps the SO(3) 
group twice. 

If the vector w does not depend on r , then the solution of equation 
(|106|) has the form 

U(t, r ) = U(t , r ) (tos ^M(t - r )^ ■ I + sin ^M(r - r )^ ■ B^j 

where U(r , r ) = ±1 and we assumed \w\ ^ . In the case |w| = the 
matrix U(r, r ) = U(r , r ). 

The usage of the SU(2) representation of SO (3) matrices not only allows 
us to reduce the number of free parameters, but is also helpful for some 
analytical calculations. Consider, for example, the matrix 



f -/d f \ t t t\ I Im{a) Im{b) - i ■ Re(b) \ 

L = , (Jfe(o) • I - U) = ^ + ^ ^ _ /m(a) j 

It is easy to verify that this matrix commutes with the matrix U, i.e. 

L = ULU* (109) 

Comparing ([1091) and (|107D we get 

Lemma 9: Let U be a SU(2) matrix, corresponding to a given A G 
SO(3) . If A ^ I then 

H = i 1 , ( /m ( 6 )' -^ e ( & ) ; Im(a)) 
a - (Re(a)) 2 



is a unit eigenvector of the matrix A, corresponding to unit eigenvalue. 



6 Normal Forms for the Spin-Orbit Hamilto- 
nian 

In this section we will consider the Hamiltonian system associated with 
the coupled spin-orbit Poisson bracket in a neighbourhood of a stationary 
point (or periodic solution) which can be canonically transformed to be the 
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origin^. In order to understand the properties of the solutions of such systems 
it is helpful firstly to find a coordinate substitution which reduces the original 
equations to the simplest possible form. Here we describe an algorithm which 
allows us to make coordinate transformations working not with equations of 
motion, but directly with the Hamiltonian function. 

6.1 Canonical Transformations of Spin-Orbit Variables 
which Map the Origin into Itself 

We shall be working with a phase space consisting of 2n + 3 variables 
z = (x, s) = (q, p, s) = (qi, ...,?„, Pi, . . . , p n , Si, s 2 , a 3 ) 



In this subsection we shall study canonical (i.e. the Poisson bracket ( 18|) 
preserving) nonsingular (i.e. with nondegenerate Jacobian matrix) maps 

z f = Z(zi) (110) 

where the components of the real vector function Z(z) are power series in 
z without constant terms. 

6.1.1 Linear Canonical Transformations 

The map which linearizes any symplectic map in the neighborhood of a 
fixed point is always a linear symplectic map. Does this apply to our case, 
where we are dealing with maps which are not symplectic but instead are 
canonical with respect to the Poisson bracket flTED? 

Lemma 1: The nonsingular real linear transformation 

z f = Azi 

is canonical if and only if the matrix A has the form 



A = diag {A orbt , A 



; '■spin J 



8 If 2*(r) is a solution of a Hamiltonian system associated with the coupled spin-orbit 
Poisson bracket then the parallel displacement z new = z Q id — z*{t) will be a canonical 
transformation if and only if s*(r) = 0. The question of what type of simple canonical 
coordinates can be introduced for an arbitrary solution is studied in Appendix B. 
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where 

A orbt e Sp(2n), A spin e S0(3) 



In more detail the map ( |110|) can be written as 

x f = X(xi, Si) 

s f = S(xi, 

Using the condition of spin-orbit Poisson bracket preservation fl50|) it is pos- 
sible to show that always 

S(x, 0) = 

Thus we have the following: 

Lemma 2: The linearization of the canonical map ( \1 1 Q ) is a linear 
canonical map if and only if 

a? °) = 

So the answer to the above question is negative in the general case. The 
linearization of the nonsingular canonical map (in its usual meaning) is not 
bound to be a linear canonical transformation. 

Definition 1: A polynomial P(z) = P(x, s) is quasi-homogeneous 

of degree m in its arguments x, s if it satisfies the equation 

P(tx, t 2 s) = t m -P(x, s) 

for every value of t. 

The set of all quasi-homogeneous polynomials of degree m {m = 0, 1, 2, . . .) 
will be denoted by TC s (m). 

Definition 2: The quasi-linearization of the map (\110j ) is defined as 
a linear transformation 

Zf = A s Zi 
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which acts invariantly on classes Ti s {m) [] and minimizes the Euclidean norm 
of the difference 



f(n) - 1. 

oz 



Lemma 3: The quasi-linearization of the map ( \11Q ) exists, is unique, is 
a linear canonical map and is defined by the matrix 

A. = diag (f (0, 3) , § (0, 0)J (111) 

Remark: The equality ( |111| ) can be used as definition of quasi-linearization. 
In this case definition 2 becomes a statement which must be proven. 

6.1.2 Factorization Theorem 

The use of the change of coordinates in the form of the power series 
( |110| ) is not too convenient for our purpose because the proof of canonicity is 



reduced to the checking of an infinite number of conditions for the coefficients 
of the Taylor expansion of the vector function Z(z). Dragt and Finn [18|] have 



shown that any symplectic map in the neighborhood of a fixed point may 
be represented in the form of a composition of a linear symplectic map and 
a sequence of displacements along trajectories of the Hamilton system with 
Hamiltonians which are homogeneous polynomials of powers 3,4,5,.... A 
similar factorization is applicable in our case too, if instead of homogeneous 
polynomials we use Hamiltonians from the classes H, s (m). 

Definition 3: We will say that the function V(z ) £ O s (m) , if 

3 lim — V(tx, t 2 s) £ H s (m) 
t^o t m 

Factorization Theorem A: For every canonical map of the form t \110j ) 
and for every integer m > 3 functions Fk(z) £ Ti. s (k) (3 < k < m) can 
be found, such that 

Z(z) = m : A s : exp(: F 3 :) . . . exp(: F m :) z 



9 This means that for an arbitrary m = 0, 1, 2, . . . and for an arbitrary P(z ) G H s (m) 
we have that P(A s z) £ TL s (m). 
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Here A s is the matrix of the quasi-linearization of the map (\11U{ ) and the 
symbol — m ' means that the difference of the right and left parts is a function 
from the class O s (m) . 

Here the exponential Lie operators are denned as usual as a power series 
(see, for example |19| ) involving the spin-orbit Poisson bracket fll8|) , and the 



action of the operator : A s : is defined by means of the rules 

: ,1, : f(z) = /(: A s : z) 
for a smooth function /, and 

i A s z = A s z 

for the identity mapping z. 

The proof of the factorization theorem A is outlined in Appendix C. 

Applying theorem A to the canonical map Z' Y {z) and then inverting the 
factorization obtained we get 

Factorization Theorem B: For every canonical map of the form 
and for every integer m > 3 functions Gk{z) G 7i s (k) (3 < k < m) can 
be found, such that 

Z(z) = m exp(: G m :) . . .exp(: G 3 :) : A s : z 



Remark: The classical Poisson bracket of two homogeneous polynomials 
of degree m and I is a homogeneous polynomial of degree m + / — 2. The 
same property holds for functions from classes TC s (m) and TC s (l) with respect 
to the spin-orbit Poisson bracket. 

6.2 General Description of the Normal Form Algo- 
rithm 

Consider a periodic real Hamiltonian function H(t, z) of period 2tc with 
respect to r which in a neighbourhood of the point z = can be expanded 
in quasi- homogeneous polynomials in z beginning with degree 2 

H(r, z) = H 2 (t } z) + H 3 (r, z) + . . . + H k (r, z ) + . . . (112) 

The purpose of this subsection is to find the simplest possible form (normal 
form) to which the Hamiltonian function (|1 1 2|) can be reduced by means 
of a nonlinear canonical coordinate transformation of the type (|110| ) which 
depends periodically on r as a parameter. 
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6.2.1 Extension of Phase Space 

Introduce two additional variables E, v and define the Poisson bracket 
and the new Hamiltonian by the rules: 

H(v, E, z) = H(v, z) — E 



{E, v} = 1, {v, Zi} — {E, Zi} — 0, i = l,2,...,2n + 3 

According to the factorization theorem B, for any given truncation order 
m every canonical transformation of coordinates of the type ( |110| ), depend- 
ing periodically on r as a parameter, can be represented with the required 
precision in the form of an operator 

3 

J] cxp(: G k (r, z) :) : A s (r) : (113) 

k=m 

In order to extend ( [L13| ) to the canonical transformation in the extended 
phase space, we rewrite ( |113| ) in the form: 

3 

IJ exp(: G k (v, z) :) : A s (v) : (114) 

k=m 

and define the action of the operator : A s (v) : on the identity mapping 
(v, z, E) by the rules 

: A s (v) : v = v, : A s (v) : z = A s (v) z 



: A s (v) : E — E — - Aj rbt J A orbt x ■ x + - s ■ cmlg (Aj pin A spin s 
Here 

A s (v) = diag(A orbt (v), A spin (v)) 
a 1 _ dA orb i i dA S pi n 

and the 2n x 2n matrix J is the symplectic unit matrix. 
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The exponential Lie operators are denned again by their power series 
with the help of the extended Poisson bracket for which we will keep the 
same notation {*, *}. 

Now, informally speaking, the problem of finding the normal form can 
be expressed as follows: Find the simplest quasi-homogeneous polynomials 
Hk € 7~C s (k) (2 < k < m) such that the initial Hamiltonian function 7i can 
be reduced to the Hamiltonian 



m 



H = > . H h (v, z) + H >m {v, z) - E, H >m e O s (m + 1) 



k=2 



by means of a nonlinear canonical coordinate transformation, which is defined 
by the operator ( 114j) . Here informal means that we have not yet defined the 



exact sense of the word "simplest". The usual way to do this is to introduce 
the simplest form for the quasi-quadratic polynomial Hi axiomatically and 
then to define Hk as functions which satisfy the condition {Hk, H 2 — E} = 0. 

6.2.2 Linear Normalization 

First we wish to simplify the quasi-quadratic part H2 of the Hamiltonian 
Ti by means of the canonical transformation : A s :. Denote 

H 2 = - P(v) x ■ x + w(v) ■ s 

where P(v) is a real symmetric matrix of order 2n and P(v + 2ir) = P(v), 
and where the real vector w(v) satisfies w(v + 2n) = w(v). 

Applying the operator : A s (v) : to the Hamiltonian 7i we get the 
following formula for the quasi-quadratic part H 2 of the new Hamiltonian 

Hi = — Aj rbt (PA or bt + JA' orbt ) x ■ x + 



spin 



+ \ s -cml g (A~l pin (C{w)A spin - A' S} 

Below we will study the most important case for applications (to which 
we will refer later on as the orbital elliptical case) when all the eigenvalues 
of the one-turn revolution matrix of the system 

dx , , 

— = JP(v)x 
dv 
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are distinct and lie on the unit circle (the general algorithm will be published 
elsewhere). Combining the results of the previous section and the well known 
normal form theory of classical linear Hamiltonian systems (see, for example 
p0| , |2~T| ) we are able to find real matrices A spin e SO(3) and A or u £ Sp(2?7,) 
27r-periodic in v such that H 2 finally takes the form 



H2 = -y (ql + Pi) + A s • s 3 



where the values ae^ and X s are called the characteristic frequencies (or 
the linear orbital and spin tunes). 

6.2.3 Recursive Loop of Nonlinear Normalization 

Let the functions Hk and Gk be already defined for k — 3, .. .,1 — 1. 
Applying the operator exp(: Gi :) to the Hamiltonian 

H 2 + ... + H1-1 + Hi + ... — E = 

exp(: G1-1 :) ... exp(: G 3 :) : A s : H(v, E, z) 

and collecting the remainders from the class Ti, s (l), we get the equation for 
obtaining the functions Hi and Gi 

f)G 

—± + {G u H 2 } = Hi - Hi (115) 
ov 

This is the so-called homology equation. We say that the quasi-homogeneous 
polynomial Gi takes the quasi-homogeneous polynomial Hi into the quasi- 
homogeneous polynomial Hi if (|115|) holds. 



In order to solve the homology equation ( |115| ) we introduce complex co- 
ordinates w = (ff, £, u) related to the old coordinates 
2 = {q-> Vi by rneans of a standard linear transformation z = Qw : 

1 + i , . N 1 — % , 
Qk = —^-{Vk + 4fc), Pk = TT-vlk ~ W 



k = 1, . . . , n 



1 + i 1 — % 
si = —z—(ui + u 2 ), s 2 = — («i - u 2 ) 
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S3 = U3 

It is easy to verify that the matrix Q has the following properties: 

QQ* = J, Q 2 = -Q* 

where the asterisk "*" indicates complex conjugation of a matrix. 

This is a Poisson transformation, but it is not canonical (although it is 
symplectic with respect to orbital variables). The nonzero Poisson brackets 
of the new basis functions are given by the equalities: 

{Vk, = 1, u 2 } = «3, {ui, u 3 } = —iui, {u 2 , u 3 } = iu 2 

In the new coordinates the equation (|115|) has the form 

^ + {gu M = h-hi (H6) 
where we have used the notation 

h 2 {w) = H 2 (Qw) = i&kTlkik + \sU 3 , hi(v, w) = Hi(v, Qw) 

g t (v, w) = Gi(v, Qw), h t (v, w) = H t (v, Qw) 

Expanding the coefficients of the monomials in the polynomials gi , hi and 
hi in Fourier series with respect to v ™ 

9i = 9i fUL exp(iA^f) rf^u 1 
< hi = hf IJL exp(iNv)7] I ^ J u L (117) 

~h { = h^ IJL ey^(iNv)ri I i J u L 

(here the summation is made over all integers iV and over all nonnegative 
integer vectors 7, J G R n , L e R 3 satisfying the condition |/| + | J|+2|L| = Z[^]) 
and substituting them in the equation ( |116| ) we obtain a system of equations 
for the coefficients: 

i ■ [N + (I - J) ■ ae + (L 2 - L x ) ■ \ s ] ■ g? IJL = hf IJL - ~hf IJL (118) 

10 We assume that this Fourier series converges absolutely. 

1 1 Here the function | * | for the integer vector k G R m is defined as \k\ = \ki | + . . . + 1 k m \ . 
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Definition 4: The characteristic frequencies \ s satisfy a 

resonance relation of order K if there exist integers ki not all equal to 
zero such that 

kieei + ... + k n as n + k n+ i\ s = (mod 1) (119) 

| fci | + . . . + | k n | +2-| k n+1 I = K (120) 



The number of the linearly independent integer solutions k\, . . . , k n+ i of 
equation ( |119| ) is an important characteristic of the frequencies aei, . . . , ae n , X s , 



and we will call it the multiplicity of the resonance. 

Note that our definition of the order of a resonance (multiplier 2 in front 
of k n+ i in ( |12C1| )) is different from the usual one, and corresponds to the 
definition of quasi- homogeneity. 

Definition 5: A quasi-homogeneous polynomial 

hi(v,w) = /^ /JL exp(iiVt;)7/eV 

is said to be a complex normal form if 

| iV + (J - J) • ae+ (L 2 - Lx) ■ X s | • | h? IJL | = 

that is, the normal form contains only resonant terms. 
In each quasi-homogeneous polynomial 

h t = hf IJL exp[iNv)n I i J u L (121) 

we isolate the resonant part h' t , which contains all those and only those 



terms of ( 121 ) whose indices satisfy 

N + (I - J) ■ ae + (L 2 - L x ) • A s = 

A quasi-homogeneous polynomial h" = hi — h[ is called the nonresonant 
part of the quasi- homogeneous polynomial hi. Obviously, hi = h[ + h'[. We 
also introduce the resonant H[ and nonresonant H" parts of a quasi- 
homogeneous polynomial Hi(v, z) as 

H&v, z) = h'i(v, Q*z), H'/(v, z) = h?(v, Q*z) 
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where h[(v, w) and h'{(v, w ) are the resonant and nonresonant parts of the 
quasi- homogeneous polynomial hi(v, w) = Hi(v, Qw). 

According to the definition 5, if the quasi-homogeneous polynomial hi is 
a complex normal form, then h[ = hi . 

Definition 6: A quasi-homogeneous polynomial 

hi(v,w) = hf IJL e^{iNv)n I i J u L 
is said to be a nonresonant complex normal form if 

|| N | + | I - J | + | L 2 - Li || ■ | hf IJL | = 

that is, the normal form contains only trivial resonant terms. 

Definition 7: A quasi- homogeneous polynomial Hi(v,z) is said to be a 
(nonresonant) normal form if after the standard transformation z — Qw 
we obtain the quasi- homogeneous polynomial hi(v, w) = Hi(v, Qw) which is 
a (nonresonant) complex normal form. 

It is easy to see (comparing definitions 5 and 7, and equation ( |118| )) 
that for every given quasi-homogeneous polynomial Hi there is the quasi- 
homogeneous polynomial Gi which takes Hi into the normal form Hi. Here 
the nonresonant part G" and the normal form H\ are uniquely determined, 
and the resonant part G\ can be specified arbitrarily. 

Among the quasi- homogeneous polynomials Gi which take Hi into the 
normal form Hi we select the one for which G\ = and we call it basic. 

Remembering that the initial Hamiltonian function 7i is real (the values 
of the function are real for real values of arguments) and that after a real 
canonical transformation : A s (v) : we again have the real Hamiltonian func- 
tion : A s {v) : Tl, we wish to find conditions on the functions Gi that will 
guarantee the reality of the transformation from : A s (v) : TC to normal form 
H. Denote by h(v, w) the quasi-homogeneous polynomial in which all the 
coefficients of the powers of w are complex conjugates (as functions of v) of 
the corresponding coefficients in the quasi- homogeneous polynomial h(v,w). 

Lemma 4: Suppose that in the homology equation (\Ll3i) Hi and Hi are 
real functions. Then the normal form Hi will be a real function and the 
necessary and sufficient conditions that the quasi-homogeneous polynomial 
Gi (which takes Hi into the normal form Hi ) will be a real function are that 
its resonant part G[(v, z) satisfies the condition 

gi(v, w) = g[{v, -Q*w) where g[(v, w ) = G\(v, Qw ) (122) 
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In particular, the basic solution of the homology equation Gi (G[ = 0) which 
takes Hi into the normal form Hi will be a real function. 



6.2.4 Normalization Theorem 

Now we are ready to summarize the results of this section as follows 
Normalization Theorem: In the orbital elliptical case, for a real initial 
Hamiltonian of the form 

oo 

H(v, E, z) = J2 H x(v, z) - E, H k E H s {k) 

k=2 

and for every integer m > 3 there exists a real canonical transformation 
depending periodically on v and defined by an operator 

3 

Y[ exp(: G k (v, z) :) : A s (v) : 

k=m 

under which the initial Hamiltonian becomes a real Hamiltonian 

ap m 
H(v, E,z) = -y (q 2 k +pt) + X s ■ s 3 + Y,H k (v, z) + H >m (v, z) — E 



k=3 



such that H >m G O s {m + 1) and all quasi-homogeneous polynomials 
H k G 7i s (k) are normal forms. Here the resonant parts of the quasi- 
homogeneous polynomials Gk may be chosen arbitrarily satisfying condition 
(\122{); then the nonresonant parts of Gk and Hk are uniquely defined. 



7 General Properties of Systems whose Hamil- 
tonians are Normal Forms 

A Hamiltonian system whose Hamiltonian function is a normal form has, 
as a rule, an abundance of symmetries, and this enables the order to be 
lowered. This is the advantage of a normal form over an arbitrary Hamilto- 
nian function. We will not discuss the technical details of the procedure of 
lowering the order here, but concentrate on the question of finding invariant 
functions (conservations laws). For convenience we will use the complex co- 
ordinates w, but it is helpful to remember that in the original coordinates z 
all functions that we use will be functions with real values. 
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7.1 Normal Forms and Invariant Functions 

Writing a normal form Ti in complex coordinates w we have 

oo 

h = h(v, E, w) = H(v, E, Qw) = iee k r] k £, k + X s u 3 + ^2h k - E (123) 

k=3 

where 

h k = h% IJL ex.v(iNv)r} I £ J u L 

and the summation is made over all integers N and over all nonnegative 
integer vectors /, J, L satisfying the conditions 

|/| + |J| + 2\L\ = k 

N + (/ - J) • ae + (L 2 - L x ) • \ s = 

We denote by £ the linear subspace of R n+l that is the linear hull of 
all integer solutions of equation ( |119| ). The dimension of C is equal to the 
multiplicity of the resonance. Let a vector m e R n+1 be orthogonal to C, i.e. 
fa _L L. Then the function 

Frn = im k r]kik + m n+1 u 3 (124) 

satisfies the condition {F^, h} = and is hence an invariant. So we have 
n + 1 — dim£ functionally independent invariants of the type ( |12 
Besides that an arbitrary normal form commutes with the function 

F = h - {i^hVhik + A s « 3 - E) (125) 

and with the Casimir function 

s i + s l + s i = 2iuiu 2 + u\ (126) 

Remark: And, of course, the Hamiltonian ( |123|) is always a constant of 
motion (in extended phase space). 
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7.2 Action-Angle Type Coordinates 

Instead of the standard cartesian coordinates 

r)i,...,r) n , £i,...,f„, ui,u 2 ,u 3 

for normal form analysis and construction we can use action-angle type 
canonical polar coordinates 

Pi, ■ ■ ■ ,Pn+2, <fi,...,<f n +i (127) 

introduced with the help of the formulae 

1+i / — . 1 + i I — 

Vi = r v2ftexp(i^), & = — — J2p l exp(-up l ) 



I = 1, . . . , n 



ui = ——^2p n+1 exp(-iip n+1 ), u 2 = ——^2p n+1 exjp(i(p n+1 ) 

U 3 = Pn+2 

Note that this change of variables is not defined on the set 

Vl ■ ■ ■ ■ ■ Vn ■ £l • • • • ■ Cn • U X ■ U 2 = 

and, as in the case of the standard linear transformation Q, it is a Poisson 
transformation, but it is not canonical (and is not even symplectic with 
respect to the orbital variables^). 

The nonzero Poisson brackets of the new basis functions are given by the 
equalities 

Wh Pi} = 1, / = 1, . . . ,n 

{Pn+l, Pn+l} = Pn+2, Wn+1, Pn+2} = 1 



1.2 



The symplccticity with respect to the orbital variables can be easily achieved, if nec- 



essary, by interchanging the places of the orbital actions and angles in (127) 
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and the Hamilton equations with the Hamiltonian function H(p, (p) take the 
form 

dpi dH dip i dH 

dt d<fi ' dt dpi 



I = 1, ... ,n 



dp n+1 dH dip n+1 dH dH 

Pn+2 ■ Tj > T7 — — — Pn+2 ' ~7\ + 



dt d(p n+1 dt dp n+1 dp n+2 

dpn+2 dH 
dt d(p n+ i 

The direct transition from the original variables z to the polar coordinates 
introduced above has the form 

qi = y / 2/^sin(v?i), Pi = yfepi coa(,<Pi) 
I = 1, . . . , n 

si = \]lp n+l cos(v9 n+ i), s 2 = y/2p n+1 sm(ip n+1 ) 

s 3 = Pn+2 

7.3 Normal Form in Nonresonant Case 

In the nonresonant case (dim£ = 0) the normal form ( |123| ) contains 
only those monomials for which 

I — J, Li — L 2 , N = 

and hence we can rewrite it as 

h = i&kVktk + X.u 3 + hi Ll L 3 (vO 1 (uiu 2 ) Ll u% 3 - E (128) 
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where the summation is made over all nonnegative integer vectors / and over 
all nonnegative integers L\, L 3 satisfying the condition 



2\I\ + AL X + 2L 3 > 3 
According to ( |124| ) the functions 

irjiZt, irj n Cn, u 3 (129) 

are constants of motion. Taking into account that h is independent of v, 
and subtracting from the Casimir function ( |1 26| ) the invariant u\ we get two 
additional invariants 

E, iu\Ui (130) 

Altogether (|129|) and (|130|) give us n+3 independent integrals of motion. This 
allows us to consider the system with the Hamiltonian (|128| ) as integrable. 
Moreover it is easy to find the general solution by quadrature. Let us do this 
in the original real variables. 

In the variables z the Hamiltonian (|128|) can be written as 

H = H(J u ...,J n} h, h) - E (131) 

where 

J k = l,(qt+P 2 k), h = ^(si + sj), h = s 3 

are the functions ( |129| ), (|130|) expressed in variables z, and the equations of 
motion take the form 

d( lk dp k 

— = UJ k p k , -r- = -UJkQk (132) 

dr dr 



dsi ds 2 ds 3 

— — = -il s 2 , — = ilsi, — = (133) 
dr dr dr 



Here we have used the notations 



dH 

u k ( Jj, . . . , J n , ii, J 2 ) = t— , (134) 
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Q(J lt 



dH 



dH 



(135) 



and we have omitted the trivial equations for the variables v and E. 

Fixing for simplicity r = and taking into account that the functions 
(134) and (|135|) are constants of motion, we get the solution of the system 
([TBI ), ([1331) in the form 



q k (r) = cos(w°r) q k (0) + sin(u;°r) p k (0) 
Pk{r) = -sin(^r) q k (0) + cos(^°t) p fc (0) 
Sl ( r ) = cos(fi r) si(0) - sin(fi°r) s 2 (0) 
s 2 (r) = sin(Q°r) si(0) + cos(ftV) s 2 (0) 



«3 r 



*b(0) 



where c<j° and f2° are the values of the functions ( |134| ) and ( |135| ) calculated 
at the initial time r = . 

Obviously, these formulae give us the solution of the triangular system 
too, if we redefine u>k and Q in (|134| ), ( |135| ) as 



(Jl, ■ ■ ■ , Jn) 



dH 



dJu 



f2 (J\ j • • • ; Jn 



h=h=o 



dH 



h=h=o 



The values iOk and Q are called the (nonlinear) orbital and spin tunes 

respectively. 



7.4 Normal Form in the Case of a Single Isolated Res- 
onance 

We shall say that the Hamiltonian function Ti is a generalized single 
resonance normal formf^] if there exist integers 



o 



N.+i> 



N° 



13 In order to keep the size of this paper within reasonable bounds we shall not discuss 
the formal aspects of constructing of the real canonical transformation which will bring 
the initial Hamiltonian to a single resonance normal form because this procedure is almost 
the same as that described in subsection 6.2. 
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where not all kf are zero such that when written in complex coordinates w 
this function has the form 

h = iaskriktk + \ s u 3 + h NIJL exp(iNv)rj I C J u L - E (136) 

and so that in the expansion (|136| ) there are only those terms whose indices 
for some integer m = m(N, I, J, L) satisfy 

Ii - Ji = m • k°, I = 1, . . . ,n 

L 2 - Lx = m-k° n+1 , N = m-N°, \I\ + | J\ + 2\L\ > 3 
Thus, the normal form ( |136| ) contains only those monomials for which 
N + (J — J) • ae + (L 2 - Li) • A s = 

= m ■ (iV° + • aei + . . . + k° ■ ae n + k° n+1 ■ A.) 
and we do not require the condition 

N° + k° ■ ae! + . . . + k° n ■ ae n + k° n+l ■ X s = 
to be satisfied. 

Let a vector m e i? n+1 be orthogonal to the vector k° . Then the function 

Fm = im k r] k £ k + m n+1 u 3 (137) 

is an invariant. So, together with the Casimir function ( |126|) and the Hamil- 
tonian ( |136| ), we have n + 2 functionally independent integrals of motion, and 
in the remainder of this subsection we shall find an additional invariant that 
will allow us to consider the single resonance problem as integrable. 

7.4.1 Single Resonance Between Orbital Tunes 

If we have a single resonance between orbital tunes (i.e. k^ +l = 0) then 
the required additional invariant is given by an expression 

F = ieklr] k ^ k + (h + E - iee k T] k ^ k - A s w 3 ) (138) 
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where the quantity 



_ as k k° k + N° 
II k° || 2 

is called the distance from orbital resonance. If the integral (|138[) is 



independent of the integrals of the form (|137|) and of the Casimir function 
(|126f) , the single resonance normal form ( |136|) is completely integrable. If 
not, then h + E is a series in irji^x, . . . , iT] n ^ n , iuiu 2 , u% only and the addi- 
tional integral will be F = «A;°%Cfc- Hence, we are always in the situation of 
complete integrability. 

The triangular truncated equations of spin motion corresponding to the 
Hamiltonian ( |136j ) 



dm du 2 .„ du 3 dh 

-i\lui, —— = t\lu 2} ~3 — = 0, \l 



dr dr dr du% 



(139) 

Ul=U2= u$= 



can be easily integrated by quadrature if we know the solution of the orbital 
part of equations of motion. In the original variables the system ( |139| ) will 
take the form 



dsi ds 2 ds 3 dH 

— = -\ls 2 , -3— = S2si, — = 0, II = — 
dr dr dr 0S3 



S1=S2=S3=0 



and this allows us to draw a number of conclusions. For example, we can 
state that in principle, it is possible to organize slow extraction of a polarized 
proton beam using a third-integer orbital resonance without significant loss 
of polarization. 

Remark: Note that the case considered here includes, in our opinion, 
the correct treatment of the situation when the synchrotron frequency is very 
small compared to the other frequencies in the system (as, for example, in 
the case of HERA-p ring |IT||). 



7.4.2 Single Spin-Orbit Resonance 

If fc° +1 7^ then as an additional invariant we can take the function 

F = eu 3 + (h + E - iee k r]k£k ~ A s m 3 ) (140) 
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where the quantity 



ae fc fc° + X s k° n+1 + N° 

is called the distance from spin-orbit resonance. If the integral ( [L40|) is 
independent of the integrals of the form ( |137| ) and of the Casimir function 
(|126|) , the single resonance normal form ( |136|) is completely integrable. If not, 
then h + E is a series in irj^i, . . . ,ir) n £ n ,iuiU2,u 3 only, and the additional 
integral will be F = u 3 . Hence, again we are in the situation of complete 
integrability. 

In more detail the invariant ( |140| ) can be written as 



F(v, ff, £ u) = eu 3 + h NIJL exp(iNv)r] 1 ^ J u 



JtJ„.L 



and we see that in general it does not satisfy the condition F(v,rf,^,0) = 
0, which is very important for application to the triangular system. This 
situation can be easily changed if we choose n invariants of the type (|137| ) in 
the form 



U 3 , 



v„ 



U 3 



and instead of F use the invariant 



(141) 



F- J2h 

\i\>i 



OIIO 



V 1 



So the triangular system will have the spin dependent integral of motion 



dF 



Ui + 



dF 



u=0 



dm 



u 2 + 



dF 



u=0 



dm 



■ u 3 



u=0 



and hence the solution for spin can be obtained by quadrature following, for 
example, the procedure described in Appendix D (a detailed analysis will be 
published in a separate paper). 



Remark: The set of invariants (|141|) show us that a proton beam cannot 
be split into two polarized parts with different amplitudes using the classical 
Stern-Gerlach effect at a single spin-orbit resonance as has been suggested in 



||34|| . For related comments see [^, [35 
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8 Methods of Numerical Integration 



Today symplectic tracking methods for orbital motion (methods which 
conserve the classical Poisson bracket^) are common tools in accelerator 
physics and we believe that such methods, which conserve the Poisson bracket 
fll8| ) exactly or with high accuracy, are the most suitable ones for numerical 
simulation of the equations (pTD-( fffi ). For simplicity we consider the case 
of autonomous Hamiltonian functions. In the nonautonomous case we can 
repeat all the steps, introducing two additional canonical variables to obtain 
an autonomous Hamiltonian system in a higher dimensional phase space. 

Those methods, of course, are mainly of theoretical interest, but in appli- 
cation to the triangular system (practical interest) the Hamiltonian approach 
allows us to reduce the initial problem of numerical integration of the system 
d5)-(0) to that of symplectic integration of the equations of orbital motion 
i) only. 

8.1 Canonical Numerical Methods for Hamiltonian Equa- 
tions 

Let K be a compact subset of the manifold M, and the autonomous 
Hamiltonian function H(z) be zero onto the set M \ N (here iV is some 
compact subset of M containing K). Then the solutions of the Hamiltonian 
system ([ll]) generate the one-parameter group of the canonical transforma- 
tions of phase space 

T(r) : M -> M, T(0) = / 

Consider a one-parameter family of canonical transformations 

T m (r) : M - M 

defined for | r | < r (0 < r ). 

We will say that T m (r) approximates T(r) in a neighbourhood of the 
identity mapping with the order m (with respect to K), if for all z E K 

II T(t)z- T m {r)z\\ = 0(r m+1 ) 

14 Numerical integration schemes for differential equations which preserve other qualita- 
tive properties (like autonomous Hamiltonian or, more generally, first integral conserving 
algorithms, and so forth) are discussed in |22], and more recent references may be found 
in || 
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Here || • || is the norm in Euclidean space, containing the manifold M. 

A one-parameter family T m (r) is said to be locally computable if for 
each zq 6 K there is a chart U (from the atlas of the manifold M) containing 
all points 

z T = T m (r) Zq 

for | r | < r and a vector function F\j such that z T is defined uniquely by 
the equation 

Fu(T, Zq, Z t ) = 

Definition: A locally computable one-parameter family of canonical 
transformations T m (r) is said to be a canonical integrator of order m of 
the group T(r), if T m (r) approximates T(r) with the order m. 



8.1.1 Recursively Generated High Order Canonical Integrators 



Let us assume that by some method we find a canonical integrator of 
order m of the group T , which has in each of the charts U the formal 
representationF] 



T m (r) = exp(: -r H + r m+1 R m (r, U) :) + 0(r 



jr+l\ 



(142) 



where m < r. 

The purpose of the present section is to demonstrate a general approach 
for constructing higher order integrators starting from integrator ( |142| ) of 
order m. 



Following 1 24, 25] we consider the mapping 



P(t) = II TSlfaT) 



(143) 



where the exponents qj are nonzero integers and the coefficients pj are 
real numbers. Substituting the representation ( 142 ) in (|143|) and combining 



15 The representation (142) holds if, for example, we assume the constancy of the rank 



of the Poisson bracket, the possibility of introducing Darboux coordinates globally in each 
of the charts U and some smoothness properties of the integrator considered. 
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the product of exponential operators into a single Lie exponent using the 
Campbell-Baker-Hausdorff formula we have 

P( r ) = 0(r r+1 ) + 

+ exp(: -r(q jPj )H + T m+ \q jP f +l )R m ^U) + 0(r m+2 ) :) (144) 
If qj , pj satisfy the conditions 

q rPj = 1, q y pf +1 = (145) 



and m + 1 > r then from ( |144[ ) it follows that the mapping ( |143| ) is a 



canonical integrator of order m + 1 of the group T(r) . 

Besides that, if m is an even number, m + 2 > r, T m (r) is time 
reversible, so that 

T m (r)T m (-r) = / 

and 

Qn+l-j — qj, Pn+l-j — Pj 

then ( |143| ) gives us a time reversible canonical integrator of order m + 2. 

Obviously, the scheme (|143|) may be recursively applied any number of 
times to obtain a canonical integrator of required order I (of course, with 
the assumption that this order I is not bigger than r). 

Example : Let 

( i \ —l i 
Tm i ^ ~J J I > "m ^ ~J J Q"m 

where nonzero integers r m , l m (r m ^ =W m ) f° r odd m satisfy the condition 

r / < 

' m <"in ^ w 

Then 

r . n m+l I / Z,m+1 _ n 
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so that ( P-45|) is satisfied for n = 2 and hence the integrator of order m + 1 
can be chosen in the form 



T m+1 (r) = T^(a m r)Tt(b m r) 

If additionally m and \r m \ are even numbers, and T m {r) is time reversible, 
then the integrator of order m + 2 can be taken as 

^m+2(i") = Tm 2 (a m r) T^ibmr) Tn? (a m r) 

Remark : Note that the recursive scheme considered is slightly more 
general than that used in P^L BSj. Here it is not necessary that the integers 



Qj m (EH) are positive and the initial integrator is not assumed to be time 
reversible. 

8.1.2 Low Order Canonical Integrators for Equations of Classical 
Spin-Orbit Motion 

In order to be able to start the recursive procedure described above, we 
now turn to the problem of constructing low order canonical integrators for a 
Hamiltonian system. We will give a few examples, and other helpful recipes 
(and references) may be found in [[22], [27], |2j| . 

Example 1: If J the structure matrix of Poisson bracket is constant, 
then for any constant matrix A , satisfying the condition J = AJ + JA T , 
the mapping 

z T = z + r (p{A z T + (I - A) z ) 

where 

0{z) = J grad s H(z) 

is a canonical integrator of first order. In the special case when A — |J , this 
integrator is a second order approximation, is time reversible, and is known 



as the mid-point rule [22 



Example 2: The property that the integrator is canonical is preserved 
under Poisson transformations. This simple remark combined with example 
1 allows us to construct first or second order canonical integrators for the 
equations of classical spin-orbit motion, because we know the explicit form 
of the transition formulae to Darboux coordinates. Note that although this 
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approach is applicable for arbitrary dependence of the Hamiltonian function 
on orbit and spin variables, its computer realization requires the use of more 
than one (at least two) local coordinate systems. 

Example 3: Consider a Hamiltonian function which depends on spin 
variables only in the linear combination % ■ s 



H = H(x, %■ s) 



Without loss of generality we can assume that 



1. Supplement the 



vector i with two unit vectors j and k satisfying the condition 



j x k 



so that the triplet i, j, k forms an orthogonal coordinate system. 
Then the mapping given by the system of equations 





(!) 


T 


S T = 


I 


( 




\k ) 





1 

cos(ar) — sin(ar) 
sin(a;r) cos(ar) 





( n 




) 


i 


so 




\k ) 





- x T + x -f _ 

X T = Xq + Ttp\ , I ■ S 



is a time reversible canonical integrator of second order defined in the initial 
variables (global ones in 9-dimensional phase space). Here J is the symplectic 
unit and 



(p (x, i ■ s ) = J gva,d £ H (x, i ■ s^j 
u (x, i • s^j = i ■ grad s H (x, i • s^j 



/ x T + x - ^ 
a = to [ — , i ■ so 



Example 4: If the Hamiltonian function H(x, s ) can be represented in 
the form of the sum 

H(x, s ) = Hi (x , h ■ s ) + . . . + H m (x, i m ■ s ) 
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then the mapping 

where A l H is the integrator of example 3, gives us a time reversible canonical 
integrator of second order. It solves (at least in theory) the problem of 
canonical integration for Hamiltonians linearly dependent on spin. 

8.2 Approximately Canonical Numerical Methods for 
Short-Term Tracking (Numerical Methods for the 
Triangular System) 

The length of the vector s is commensurate with Planck's constant k. 
If we perform the renormalization 

h _ 

Sold > 2 ' ^newy 

then the equations of spin-orbit motion corresponding to the Hamiltonian 
function H(x, s) become 

dx 

— = J grad g H orbt (x) + 0(h) 



= 1 



d,S r 



dr 



W(x) x s new + 0(h) 



where 

H r b t(x) = h(x,6) and W(x) = (grad g H(x, s))| s - =5 

Consequently both the effect of spin on the orbit motion and the induced 
nonlinear influence of spin on itself, are very small and we will neglect them, 
if the integration interval T satisfies the condition 

T « h- 1 



This means that with high accuracy we can reduce the initial problem of nu- 
merical integration of the equations of spin-orbit motion to that of symplectic 
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integration of orbital motion equations only. We will now demonstrate this 
in a few steps. 

First, for any vector U{x) we can approximate the action of the mapping 



exp (- : U(x ) ■ s + O (Jsf 



on phase space variables with required precision by means of the explicit 
formulae 



X j X{ 



(146) 



s f = U + 



sin 



(\U\) 1 - cosfl U \) „ 
' lJ -A(x l ) + J A 2 (x t )\ s t (147) 



U 



where 



A = 



U 



( o -u 3 u 2 \ 

U 3 -Z7i 

V c/i o y 



Second, for any given positive integer k and using the Campbell-Baker- 
Hausdorff formula it is possible to find vectors 

W\x, r), W 2 (x, r), W 3 (x, r), W\x, t) 

so that the following decomposition formulae obtain 

exp (- :tH:) = exp (- : r (H orht + W -s + O (jsf 



exp (- : ^# orbt exp (- : r (W 1 ■ s + O (\s\ 2 )) :) exp (- : T -H orht 
= * exp (- : I (w> 2 • s + O (\s\ 2 )) :) exp (- : rtf or6t :) • 
■e W (-: T -(w 2 .s + 0(\s\ 2 ))-) = k 



71 



exp (- : rH orbt :) exp (- : r (w 3 ■ s + O (\s\ 



- k exp (- : r (w 4 ■ s + O (|sf )) :) exp (- : tH, 



orbt 



(148) 



Here r is the size of the integration step and = k indicates that the differences 
between the right and left side have at least the order 0(\ r | ). 

Finally, let us select a decomposition formula from ( 148j) and hence one of 
the vectors W % (x), i = 1,2, 3, 4. If we use the combination of some symplectic 
integration method of order k for orbital motion and the formulas ( |146|) , ( |147|) 
for the mapping 

exp (: A t) -s + 0(\s\ 2 )) :) 

(A = r/2 for i = 2, or A = r for i — 1, 3, 4) we will obtain a numerical method 
which has the order k, is symplectic for the orbital motion and automatically 
maintains the length of a spin vector equal to its initial value. 
Example: 



2 2 

W 1 = W + T — ■ {H orhu {H orbt , W}} T ~ 



12 



{H orbu w}xW 



w 2 = w + 



24 



[H orbt , W) x W] - T — ■ {H orbu {H orbt , W}} 



for k = 4 and 

W 3 = W + T - ■ {H orbt , W) , 



W — w — 



T 



{Horbt, W} 



for k = 2 . 

Here the notation {H orbt , U\ means the vector with components: 
{{Horbt, Ui} , {H orbt , U2} , {H orbt , U3}) 



9 Additional Transformations of the Spin-Orbit 
Hamiltonian 



One of the most important problems in the design of accelerators is the 
problem of providing electric and magnetic fields which can hold a charged 
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particle beam in a sufficiently small neighbourhood of some geometrical line 
which we already introduced in subsection 4.3 under the name of closed 
design orbitQ. One of the commonly used approaches to find a solution is to 
create electric and magnetic fields for which the given geometrical line will 
be the projection on R 3 of the six-dimensional trajectory of orbital motion 
in these fields, such that this trajectory has to be stable with respect to 
small perturbations of the initial conditions and of the values of electric and 
magnetic fields. We will now call this trajectory the reference particle, 
and all transformations of the spin-orbit Hamiltonian in this section will be 
more or less connected with this concept. 



9.1 Closed Design Orbit and Reference Particle 

Since in the final form of the equations we would like to keep the possibility 
to treat effects like misalignments of different electromagnetic elements of 
the accelerator, fluctuations in the values of electric and magnetic fields, 
influence of space charge, etc., we will assume that the reference particle 
is not a trajectory resulting from the Hamiltonian function (jS9|), but is a 
solution of the system with Hamiltonian function 

H = H or bt + • s = — aexTTg + eeyTt^ — (1 + hx + ay)iif — ~A Z (149) 

where 

6 ~ 6 ~ 

= P% ~ gAtf' Kg = Py — 



7T f = ^ - 2 m c - ^ - TTgJ 

and the components of the vector potential Ag , Ag , A z , and the scalar 
potential $ , generally speaking, are not equal to the corresponding values 

16 In fact in this paper we never use the condition for the closed design orbit to be a 
closed curve, and so the equations derived can be used in applications to linear accel- 
erators, cyclotrons and etc. Note also that in order to be more consistent in a coupled 
spin-orbit formalism the closed design orbit should be considered as a curve in the six- 
dimensional space of three orbital and three spin coordinates, but we do it in fact by 
assuming that for the spin part this curve satisfies the equation s — for all the times. 
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in the Hamiltonian function (jS^)Q. 

Now we must discuss the conditions under which the functions 

x (z), P°(z), y (z), P°(z), E (z), t (z) (150) 

will be the solution of the equations of motion corresponding to the Hamil- 
tonian (|149|) , and the conditions under which the projection of this solution 
on R 3 will coincide with the closed design orbit. 

Since x and y were introduced in subsection 4.3 as the transverse devia- 
tions from the closed design orbit it is necessary that 

x (z) = and yo(z) = (151) 

Substituting ( |I50| ) and ( |151| ) in the equations 

dx dH or bt (1 + hx + ay) _, 

dz oP x 71 f 

dy dH orbt (1 + hx + ay) _ 

dz OPy IX f 

we obtain that the momenta P2 and P° have to be 

P° x (z) = -A%(z) and F»(z) = - A%{z) (152) 
c c 

where 

A%(z) = Ajf(x (z), y (z), z, t (z)) 



A %( z ) = A s (x (z), y Q (z), z, t (z)) 
Then putting ( |150| )- ([T52|) in the equations 

dE dH orbt dt ^dHorbt 

dz dt_ ' dz dE 

17 Somctimcs it is helpful to assume that even the charge and the rest mass of the 
reference particle are different from the corresponding values in the Hamiltonian function 
© too. 
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we obtain that E (z) and t (z) have to satisfy the equalities 

dE _e fdA z \° _e_ E - e$ / 
1z ~ ~c\Wj + ^% c 2 Vdt) 



(153) 



dto 1 E - e$ 

dz 



*°f 



where 

$o(z) = $(x (z), y (z), z, t (z)), 



7t% = 



'(Eo-e^y 



(154) 



— mpC 2 



'd& 
~dt 



<9$ 
~dt 



Z = Z,t=t(j 



The two remaining equations 



dP* 
dz 



dH, 



orbt 



dx 



'd_A 
dt 



dPy 

dz 



OA, 



dt 



Z=Z,t=to 



dH, 



orbt 



dy 



give us two additional conditions which have to be satisfied too 



dA 



N 



c dz 



1 c \ dx 



e E - e$ ( 

c 2 1 dx 



(155) 



Here 



c dz 



e 
c 



k dy 



e E n 



5 W 



(156) 



9a; 



9$ 



2!=^0>!'=2'0. 
Z=Z,t=to 



<9y / <9y 



Z=Z,t=to 



'dA 
dx 



dA, 



dx 



x=x ,y=y , 
Z=Z,t=t(, 



'dA, 



dA, 



dy J dy 



x=x ,y=y , 

z=z,t=to 
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So we find that the necessary and sufficient conditions for the functions ( |1 50| ) 
to be a reference particle with respect to a given closed design orbit and elec- 
tromagnetic field are the conditions ( |151| )-( |T56| ). These conditions (especially 
( | 15 5| ) and (|156|) ) looks rather complicated, but in practice they can often be 



satisfied with the help of a very simple model of the closed design orbit (a 
sequence of line segments and arcs) and a corresponding piecewise constant 
model of the electromagnetic fields. We will not discuss in this paper the 
conditions of stability of this trajectory (and different definitions of stabil- 
ity) because this is one of the major problems of accelerator physics and it 
deserves special detailed consideration. 

The reference particle is a trajectory, but for this trajectory just as for 
some real particle, we define quantities such as Lorentz factor 70, relative 
velocity /3q, and modulus of the velocity and kinetic momentum vectors, 
which we will denote by vq and n respectively: 



7o = — , A) = \/l 5> v o = Poc 

TJIqC 




r h F J {Eo-e<Z> ) 2 2 2 _ 
7T = m 7oW = m 7 A)C = m 7ocW 1 - = y ~ 2 mfic 2 = n f 

In these new notations the equations ( |153| ) and ( 154j) can be rewritten as 

dEo = _e fdA z \° + _e_ dto = J_ = 1_ (15?) 

dz c \ dt J (3qc \ dt ) ' dz (3qc Vq 

Note that sometimes it is more convenient to introduce the reference 
particle axiomatically just as an arbitrary set of six functions ( |150[ ) satisfying 
only the conditions 

dt 

x (z) = 0, y {z) = 0, -/ ^ 

dz 



In this case we define the value of v using the second equation in (|157|) , 
and after this j3 , 70 and tt will be defined as follows 

a v ° 1 a 

Po = — , 7o = , tto = m 7oP c 
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In order, in the following, to have a uniform treatment for both ways of 
introducing the reference particle, for the axiomatic case we will define also 
A°fi, A°g and <f> as 

o . c D o aQ . c d0 * E -m -f c 2 



= -P u = -P u $ 

9.2 Deviations of Orbital Coordinates from the Solu- 
tion for the Reference Particle 

Since the reference particle is a reflection of our understanding (or of our 
desire) of how the accelerator should behave it is natural to introduce as 
new coordinates the deviations of orbit variables from the trajectory of the 
reference particle. Let us define new coordinates (x, y, a) and new momenta 
[Pxi Py, as 



•&new ^oldi 



pnew _ pold pO _ pold e 40 

x x x x c N 



„. _ pnew _ pold pO _ pold £ 

ynew yoldi r y r y r y r y c B 

a = - (t - 1 ) , e = E - E 

The orbit Hamiltonian in new variables has the formQ 

6 

H rbt = — aex7r^ + a3?/7r^ — (1 + hx + ay) Hf — -A z + 



c 



dt dE e ( dA% dA% 

+ e — + a — — + - \x + y — ? 

az az c I az az 

where now 

= (^/V) TTfl) TTjO 



18 Note, that before this transformation E was a generalized coordinate and t was a 
generalized momentum, but now we choose a = —(t—to) as the new generalized coordinate 
and £ = E — Eq becomes a generalized momentum. 
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(22(2 i \ 2 2 \ 1 / 2 

7T f = (m Q C ( 7 -1J - 7T^ - Tig) 

AA^ = Afi — A% AA S = A s - A% A$ = $ - <f> 

e -eA$ 



7 = 7o + 2 

and the spin part keeps the same form as in subsection 4.6 with the compo- 
nent of the vector tt and 7 defined above. 

The variable a describes the difference in arrival times at the position 
z between a given particle and the reference particle, and the quantity e is 
their energy deviation. It seems to be more suitable to use another pair of 
canonical variables 

a 1 &old 

Cnew = PoCCTold an CL E new = — — 

Poc 

where the new a describes the longitudinal separation of the particle from 
the centre of the bunch. This canonical transformation does not concern any 
another variables and the new orbital part of the Hamiltonian becomes 

6 

H rbt = —sexng + eeyn^ — (1 + hx + ay) i\f — -A z + 

a dE ae dn e ( dA% dA% 
PqC dz 7r 7o dz c \ dz dz 

whereas the spin part H sp in remains the same if we take into account the 
new expression for 7 

7 = 7„ + ^-fA$) 

9.3 Scaling of the Orbital Variables 

The canonical scaling of the orbital variables which we use in this paper 



is 



pold 

t frr^T , , pnew _ x 

•^new — V "(M'Otoj 1 ; 



'vro 
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pold 

Vnew = y/^OVold, Py £W = -7= 

VTTo 



, — £oid 

O'new — y^O^old, £■ 



new 



This scaling is different from that usually used in accelerator physics (see 
Appendix E) and is applicable for both storage and acceleration regimes. 

9.4 The General Form of the Spin-Orbit Hamiltonian 
in New Coordinates up to First Order with Re- 
spect to Spin Variables 

Before writing out the final form of the spin-orbit Hamiltonian we would 
like to note that various authors use various coordinate systems for the treat- 
ment of fully coupled transverse and longitudinal motion (synchro-betatron 
motion). Our variables are closest in their physical meaning to the coordi- 



nates used in 30, 31 



In the variables introduced above the spin-orbit Hamiltonian takes the 
following final (at least for this paper) form 



H = H orbt + H_ 



spin 



Horbt = -aesf-p-J + eey(^fir-) - ( ^ + hx + ay) ( J - ~A + 



a dE e ( dA% dA° n 

+ + — - — + -— U— ^ + y — — ^ I + 

7Tq po c dz Jtiq c \ dz dz 



1 (a e x P x + y P y + a e\ dir 



H S pin = -as x + hsy - <£s z + + hx + ay) T ^l 1 W ■ s 
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where 



G(n.B)n i / i 



n x S 



m 7c I mgC 2 (l + 7) moc \ 1 + 7, 

= (^am ^B. TTf) 



^ = V^( p x &Arf , 7T 5 = ( P y AAg 



7T f = ^m2 C 2 ( 7 2 - 1) - 7T^ - 7T| 

v^oA / e \ 

7 = 7o + e A$ 

moc V V^oPoC / 

7o = n — , tto = m 7oPoc 

Remembering that the electric and magnetic fields, and the vector and scalar 
potentials are supposed to be defined in the curvilinear coordinate system 
connected with the closed design orbit (see subsection 4.6 and appendix A 
for more details) and using the conversion formulae 

x' = x/y/TrE, P' x = y^P x + (e/c)A% 

y' = y/V^o, P y = V^Py + (e/c)A% 

t> = t - a/(^(3 c), E> = E + (y^A,c) e 



z' — z 



from these variables, denoted here as z', x', P' x , y' , Py, E', t' , to our final 
canonical coordinates z, x, P x , y, P y , a, e we obtain the rule for the sub- 
stitution of the arguments of the mentioned above functions 

F(t', x>, y>, z>) -> f{u- ° ' !J - ] 
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In our final variables the relations between fields and potentials become: 
The magnetic field: 



B 



N 



+ hx + ay 



dA z dAg sey 
dy dz y/i^ 



B„ 



B 3 



'dA 



N 



y/jfo + hx + ay \ dz 



,— dA z eex 
ox Jt^ , 



where 



' dAg _ dA£ 
dx dy 



A z = 



((x/tto + hx + ay) Af + ee(xAg - y A$)) 



The electric field: 



^_ / dA s d& 



i (d$ ( <9$ <9$ 

+ + \ dx dy / 



To have the spin part of the Hamiltonian in more detailed form, let us 
introduce a vector Q by means of the equality H sp in — Cl - s and write out 
its components fl x , Q y and Q z : 



—a + (t/hq + hx + ay) 



TTq 7Tf C 



mi c 2 (1 + 7) 
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— \ G + — ) (7Tg£f - TTfSg) 



m c 



1 + 7 



^s/ = h + (v^o + hx + ay) 



TTq 7T^ C 



-{1 + 7GJ bg H 2 2n , — v h 



H (G + ; I (iTfSfi - KfiSf) 



m c 



1 + 7 



Q 2 = -ae + (0^0 + /ij; + a?/) 



7Tq 7T^ C 



-(1 + lG )B f + ±^xM + 

T m\c 2 (1 + 7) 



— IG + — — I (vr^^ - 7T 5 £#) 



m c 



1 + 7 
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A How to Transform the Equations of the 
Electromagnetic Field to Curvilinear Co- 
ordinates Associated with the Closed De- 
sign Orbit 

The purpose of this appendix is to point out the form which differential 
operators entering in the Maxwell system of equations for the electromag- 
netic field will have in the curvilinear coordinates considered (they could be 
equations for the vector and scalar potentials, or direct equations for the 

— * — * 

electric £ and magnetic B fields) 

— * 

In the coordinate system used every vector A can be uniquely represented 
in the form 

A — Ajf - N + A s - B + Ap-f 

(at least in a sufficiently small neighbourhood of closed design orbit). 
As physical Components of A we understand the values 

A^ = Afi(t, x, y, z), Ag = Ag(t, x, y, z), Af = Af(t, x, y, z) 

Introduce the vectors 

f z = (1 + hx + ay) f + ae (xB - yN^j 

N z = (l + hx + ay)N + seyf 
B z = (l + hx + ay) B - aexT 

— * 

and define A z , A x , and A y as projections of the vector A on the vectors 

— * — * — * 

T z , N z and B z respectively 

A z = A ■ T z , A x = A ■ N z , A y = A ■ B z 

Using the quantities introduced above we get the following formulae: 
The Gradient of a Scalar Function: 
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grad0 = V0 



1 + hx + ay 



dz 



T + 



jf ■ N z + ■ B z 
dx ay 



The Laplacian of a Scalar Function: 



A0 = V 2 



1 



1 + hx + ay 



«.( W .f) + |.(v#.At,) + | ; (V*-&) 



The Divergence of a Vector Field: 



divA = V-A 



1 + hx + ay 
The Curl of a Vector Field: 



dAf dA x dA y 



dz 



dx 



dy 



cm\A = V x A = 



1 + hx + ay 



'dAs 
dx 



OA 



N 



dy 



T z + 



'dA z OA ; 



dy 



dz 



N + 



+ 



dA$ dA z \~ 



dz 



dx 



The Laplacian of a Vector Field: 

The Laplacian of a vector field, if it is needed, can be expressed using the 
above formulae using the equality 



AA = V(V • A) - V x (V x A) 
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B Simple Canonical Coordinates for the Pe- 
riodic Solution 



Let the differential equations of spin-orbit motion 

dz 

£- = {z, H(r, z)} 

T-periodic in r and with possible nonlinear dependence of the Hamiltonian 
function H(t, z) on the variables s*have the T-periodic solution z*(t) 



x(t) 

satisfying the condition 



xAt) 



s{t) = s*(r) 



158) 



+ 



(159) 



Since the parallel displacement z new = z a id — 5**(r) which converts the 
solution ( |158| ) into the origin is not a canonical transformation with respect 
to the coupled spin-orbit Poisson bracket if ( |159| ) holds, in this Appendix we 
introduce other simple canonical coordinates for the periodic solution ([] 
Introduce the real skewsymmetric matrix 



C(r, z) 



-M 

OS3 



9H q 

dS3 



\ 



dH 

ds2 



dH 



dH \ 
dS2 » 



dH 

ds i 



0/ 



and write C{t) = C(r, -?*(t)) . 

It is easy to check that the unit vector 



1 



n t 



SAT 



periodically dependent on r satisfies the equation 

dn ~. . 
— = C(r)n 
dr 



(160) 
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Let us assume that we have found two unit vectors m(r) and l{r) 
T-periodic in r which supplement the vector n{r) to form an orthogonal 
coordinate system satisfying the condition 

m(r) x f(r)] ■ n(r) = 1 (161) 

Introduce new coordinates y, u with the help of the canonical transformation 

x = x*(t) + s = A(t)u 

A(t) = (m(r), f(r), n(r)) e SO(3) 
In the variables y, u the Hamiltonian function takes the form 

dx 1 / dA\ 

H new = H i d (t, f*( r ) + A(r)w) + y' J ~^T^2^' CUrltf ( ^ T dr ) 

and the periodic solution z*(r) is now expressed as follows 

y* = o, ir, = (o, o, |s*(o)|) 

Now we wish to discuss the problem of existence and possible freedom of 
choice of the vectors m(r) and Z(r) introduced above. 
Taking the derivative with respect to r in the identities 

m(r) ■ n{r) = 0, m(r) ■ m(r) = 1 

we have 

/ dm - . , ^\ . , dm 

C(r)m ■ n = and — ■ rh = 162 

\dr J dr 

Subtracting the identity Crh-rh = from the second of the equalities ( |162| ) 
we find that the vector 

dm ~, _, 
— — — C(t) m 
dr 

is orthogonal to the vectors n and rh for all values of r, and hence can be 
represented in the form 

dm ~ — » 
_ _ C (r)m = Mr) I 
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Similarly we have 



From the condition 



^ - C{r)l = Mr)m 
ar 



^L(rh-T) = 

GST V ' 

it follows that tyiir) = —i/ji(t) = f tpij) , i.e. the vectors m, / satisfy the 
system of differential equations 

^ = C(r)m-^(r)/, ^ = ^(r)Z + ?/>(r)m (163) 
ar ar 

The solution of (|163|) , written in complex notation, is 

m(r) +U(t) = exp li [ ip(rj)dr]) D(r)(m(0) + il(0)) (164) 



where -D(r) is the fundamental matrix solution of ( |160| ) so that 

^ = C(t)D, D{0) = I 

Remembering that the vectors m and I must be T-periodic in r we have 
from ([1641) 

D(T)(m{0) + iZ*(0)) = exp \ -i J ip(r])d7^j (m(0) + «T(0)) (165) 

so that 

exp ^— i y ip(i])dr] 

is an eigenvalue and m(0) + zZ(0) is the corresponding eigenvector of the 
matrix D(T). Multiplying ( |165| ) onto the vector m(0) + 2/(0) we get 

<>xp ( -i [ ip(r))dr)\ = -D(T)(m(Q) + ,7(0)) ■ (/7/(0) + /7(0)) (166) 



This arguments can be reversed to show that the differentiable unit vec- 
tors m(r) and l{r) T-periodic in r will supplement the vector n(r) to 
form an orthogonal basis satisfying the condition ( |161| ) if and only if they 
do so for r = and then satisfy the differential equations ( |163| ) for some 
continuous function ip(r) for which the equality (|166|) holds. 

Applying the transformation s = A(r) u described above to the Hamil- 
tonian 

H = w(t) ■ s, w(t + 2tt) = w{t) 

discussed in subsubsection 5.2.4 we find that in the variables u this Hamil- 
tonian becomes 

H = ip(r) ■ u 3 

and is hence a normal form for all possible functions ip(r) satisfying ( |166| ), 
which are independent of r. 



C Sketch of the Proof of the Factorization 
Theorem 

The purpose of this appendix is to prove the factorization theorem A. 
The matrix of the quasi-linearization A s is uniquely defined by the formula 
( |111| ). Introduce maps 

M 2 = : A s : , M k = M k -i exp(: F k :) , k = 3, ...,m 

Using induction in k, we will show that the functions F k G Tl s (k) can be 
chosen in such a way that the conditions 

M k 1 X(x, s) = fe+ i x + X k (x, s) 

(167) 

S(x, s) = k+ 2 s + S k +i(x, s) 

are fulfilled where X k G Ti s {k) and S k+ \ G 7i s {k + \). 

It is easy to see that the equality ( |167| ) for k = m gives the proof of the 
theorem. 



Obviously, ( |167| ) is correct at k = 2. Applying the operator 

exp(- : F t (z) :) 



to both parts of ( 167 ) for k = I — 1 we obtain that ( |167| ) will be correct for 
k = I, if the function F\ satisfies the equations 

{F h x} = X^t, {Ft, s} = S t (168) 

The structural matrix of the spin-orbit Poisson bracket has the form 

J(z) = diag(J, J B (s)) 

where the 2n x 2n matrix J is the symplectic unit and 



JJs 



( s 3 -s 2 

-S 3 Si 

\ s 2 -si 



With the help of the matrix J the equations (|168 ) can be written in the form 
of the system 



J-grad^F/ = -I w 



(169) 



J s (s)-grad s ^ = -S t (170) 
Represent the vector function X/_i in the form of a sum 



k=0 

where X^(x, s) are homogeneous polynomials of degree k in the variables 
x and the symbol [m] denotes the biggest integer which is smaller or equal 
to m. 

Lemma 1: Any solution Fi(z) G TL s (l) of the system p6'&|j -( [77^ ) is given 
by the formula 



Fi 



x 



E 

fc=0 



1 



I - 2k 



JX, 



H 

l-2k-l 



V{s 



where V(s) G TC s (l) satisfies the equation 

J s (s)-grad g V = St(0, s 



(171) 
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The proof of this and following Lemmas is left as an exercise for the 
interested reader. 

Starting from this point, we can forget about the existence of the orbital 
variables x, because the proof of Theorem 1 has been reduced to finding the 
solution of equation ( |171| ), depending only on variables s. 

Let tt s (k) be the class of homogeneous polynomials of degree k in the 
variables s. Our current task is to find the solution V(s) G H S (A;) of the 
equation 

J a (s)-gradV = f(s) (172) 

where the vector function / G H s (k), k > 1 satisfies the condition 

grad(s-/) = div(/).s (173) 

The condition ( |173j ) follows from requiring that the map ( |110| ) is canonical. 

Lemma 2: If the solution V(s) G H s (k) of the system < \172( ) exists, then 
it satisfies the equation 

k(k + l)V - \s\ 2 AV = R(s) (174) 

where 

d 2 d 2 d 2 
R(s) = s-curlf, A = M + M + ^l 



Lemma 3: If the solution of the equation ( \174 ) exists, then for odd k it 
is unique and for even k the difference between any two solutions is given by 
the formula 

V 1 {s) - V 2 {s) = c-\s\ k 

where c is an arbitrary constant. 

Lemma 4: The equation ( \1 7^ ) has a solution if and only if 

A^R(s) = (175) 



If ( |J 75p holds then the function 



V = a i\*\ 2j & j R(s) (176) 

3=0 
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where 



%j Z ^{k-2i + l){k-2i) 



satisfies the equation (\174) - 
Lemma 5: 



A j (s- curlf) = s- curl (A j f 

and therefore, 

A^l (s- curlf) = 



So in Lemmas 2 - 5 we have established that if the solution of the system 
( 172|) exists, then it is unique (up to an additive Casimir function of the 



spin-orbit Poisson bracket (|18|) c • |s| fc at even k) and is given by the formula 
(|T76|) . Substituting (|7|) in (g7g), we obtain: 

Lemma 6: The necessary and sufficient conditions for the solvability of 
the equation (\172j ) are that 

s ■ f = and divf = (177) 



The complete solution of the problem is 

Lemma 7: If f G H s (k), then conditions ( \1 77] ) and ( \1 75p are equivalent. 



Remark: From ( |167| ) it follows that the actual precision of the represen- 
tation of spin variables (functions S(x, s)) in theorems 1 and 2 is higher by 
one order (in the sense of classes O s (m)) than the precision of the represen- 
tation of orbital variables. 



D How to Integrate Linear Equations of Spin 
Motion in Quadratures if We Know their 
Partial Solution 

For general linear homogeneous systems of ordinary differential equations 
of order n knowledge of one nontrivial solution allows us to lower the order 
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of the system by one unit. So, in order to be able to find the general solution 
in quadratures we need to know n — \ linearly independent solutions. In this 
appendix we will show that due to a special symmetry, third order linear 
equations of spin motion can be integrated completely if we know only one 
partial solution. 
Let the vector 



n(r) = (ni(r), n 2 (r), n 3 (r)), 



1*001 



satisfy the system 



ds 



^ = W(t) xs = C (W(r)) s 



(178) 



Without loss of generality we will consider the r-interval on which one of 
components of the vector n(r) is not equal to zero. Let it be ^(r). If on 
the considered interval ni(r) 7^ (or n 3 (r) 7^ 0) we can begin by making a 



coordinate transformation s r 



Bs Q i d where the matrix 



B = 



/ -1 \ 

1 
V l) 




Introduce the vectors 
m(r) 



(n 2 (r), -ni(r), 0) 



l(r) = H(r) x m(r) 



(ni(r)n 3 (r), n 2 (r)n 3 (r), ng(r) - 1) 



|m(r)| = |Z(r)| = 1 



and the new variables 



u 



/ m(r) 

f( r ) ] . s d ^ A(t) ■ s 
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In the new variables the system (17?) becomes 



? = i d 4-^ + A.C(W)-A 
dr \ dr v ' 



where the components of the vector b(r) are 



u 



= C (b(r 



u 



(179) 



b, = ^-f + n-C(w)l, 
dr v ' 



, dn _ 

b 2 = — — • m 
dr 



n 



c(w 



m 



b 3 = ^-m + 1-C(W) 



dr 



Since the vector n(r) satisfies the system ( |178| ), then fei(r) = and 62 (t) = 
. After some simple manipulations we get for u = 63 



l-n|(r; 



(IVi(r)ni(r) + W 2 (r) n 2 (r)) 



(180) 



and hence the solution of the system ( |179| ) is 



it r 



/ cos ^(r,r ) - sin *(t,t ) 

sin *(r, r ) cos ^(r, r ) | • w(r ) d = K(t, r ) • w(r ) 
V 1 



where 



*(t,t ) 



So the solution of the initial system (|1 78|) can be expressed now as follows 

s(t) = A T (r)K(r,T )A(r )s(r ) 
where the matrix A(t) has the form 



A(r) 



■«2(r) 


ni(r) 






ni(r)n3(r) 


«2(T)n 3 (T) 


\/l-«i(-r) 




m(r) 


n 2 (r) 
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\ 

yi-nf(r) 
n 3 {r) J 



Remark: It is easy to see that in fact the procedure described above is 
applicable for condition weaker than n 2 (r) 7^ , n\{r) + 7i|(r) 7^ (which 
is equivalent to 1 713(7")] < 1). 

Example: The equations of spin motion generated by the model Hamil- 
tonian of the stationary single resonance problem 

H(t, s ) = e cos(fir + 0) • Si + e sin(fir + 0) ■ s 2 + A ■ s 3 (181) 

for e ytz have the nontrivial (27r/f2)-periodic solution 



rt(r) 



1 



e 2 + (A - fl) 2 



(ecos(fir + 0), £sin(fir + 0), A — Q) 



which satisfies the condition (713 (r)| < 1 for all values of r. 

The (27r/f2)-periodic matrix A(t) defined by this solution has the form 



A{r) 



( ^sin(fir + 0) -^cos(fir + 0) \ 
A - n cos(fir + 0) ^ sin(fir + 0) _M 



\e\ uj v \e\ uj 

- cos(fir + 0) - sin(fir + 



where to, calculated in accordance with (|180|), is 



u = sJ7 2 + (A - ny< 



const 



So the fundamental matrix solution of the stationary single resonance prob- 
lem 



M(r, r ) = A t (t) K(r, r ) A(t ) 
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expressed in terms of matrix elements takes the following form 
mn = asm(Qr + 0) + (^^P + (^) 2 cos(fir + 0)^ cos(fir + 0) 

m 12 = (^pP + (£) 2 cos(fir + 0)^) sin(fir + 0) - a cos(fir + 0) 
mi 3 = * cos(fir + 0) - (3) 

m 2 i = Csin(fir + 0) + (^7 + (^) 2sin ( fir + 0)) cos(fir + 0) 
™22 = (^r7 + (^) 2 sin(^r + 0)) sin(fir + 0) - Ccos(fir + 0) 
™23 = £ sin (^ + ^) - 7) 

m 3 i = ^ (-$cos(£7t + 0) - sin(cj(r - r )) sin(fir + 0)) 
= f; (sin(cj(r - r )) cos(f2r + 0) + ^sin(^r + 0)) 

. m33 = (^) 2 cos(c(r-r )) + (^) 2 

where we have used the notations 

a = cos(lo(t — t )) sin(f2r + 0) + sin(cj(r — r )) cos(f2r + 0) 

/3 = cos(cu(r — r )) cos(f2r + 0) — sin(cj(r — r )) sin(f2r + 0) 

< 7 = sin(cu(r - r )) cos(fir + 0) + cos(cu(r - r )) sin(fir + 0) 

C = sin(cj(r - r )) sin(f2r + 0) - cos(cu(r - r )) cos(f2r + 0) 

^=^0 (1 -008(^-76))) 
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E Noncanonical Scaling of the Orbital Vari- 
ables in the Case of the Simplest Storage 
Regime 

In the case of the simplest storage regime, that is 

ttq(z) = const > 0, 

it is typical in accelerators physics to keep the orbital position variables 
unchanged and to normalize the corresponding momenta to the value of the 
design kinetic momentum tt . So, following tradition, we would like to have 
the possibility to describe the orbital motion using the new variables 



•Knew •Eoldi Unew Void) @new @ old 

pnew = poldf^ pnew = pold^ ^ = 

where the superscript 'old' indicates variables of subsection 9.2. 



;i82) 



Unfortunately the coordinate transformation ( |182| ) is not symplectic and 
even if for the study of the orbital motion alone we can treat this transforma- 
tion as canonical using a rescaled Hamiltonian (see remark 2 in subsection 
4.1), we cannot proceed in the same manner for the case of fully coupled 
equations of spin-orbit motion. So we will restrict ourselves to presenting 
the Hamiltonian which will give us the correct form of the triangular system 
in new variables. This will not be applicable to the study of the effect of 
the spin on the orbit motion if we admit the complete equations of motion 



In the variables introduced above the spin-orbit Hamiltonian becomes 

Horbt Hspin 



Horbt = — — — = — sexTTg + dsyiTfi — (1 + hx + ay) tta 



e , a dE e ( dA% dA% 

— a z + e + —^-r 1 + — \ x -r L + y~r 

7T C 7T p C dZ 7T C I dz dz 
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tt i (1 + hx + ay) ? _ 
-"spin = + — aes z H 3 W ■ s 



Tlr, 



■ i . ... vr^G (vf • B ) n ^ 



where 

# = — — I (1 + 7 G)# 



7T C 



mQC 2 (l + 7) m c \ 1+7 



n x £ 



7T 



A' 



7T 



P. - 



AA 



7T C 



7Ti 



7T £ 



7T C 



7Tj! 



/ 7 2 _ 1 



^2 v 2 

— — 7T - — 7T - 

2„,2 "at "fi 



Phi 



7 = 7o + 7o/9q ^ - 



A$ 



For such scaling the conversion formulae between variables in the curvilin- 
ear coordinate system connected with the closed orbit z', x', P' x , y', Py, E', t' 
and our final coordinates z, x, P x , y, P y , a, e have the form 



x' 



x, 



P' 

x 



y = y, 

f = t - <7/(A)C), 



n P x + (e/c)A% 
P' y = n P y + (e/c)A% 
E' = E + (it {3 c)e 



z = z 

The rule for the substitution of the arguments reads now as 



F(t', x', y', z') -> F 1 1 - — , x, y, z 
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The expressions for the projections of the magnetic field and for A z keep 
the same form as in subsection 4.6, and the connection between the electric 
field and the potentials becomes 

N dx da 

<9$ „ dA 



£r = — r- + A) 



B 



dy da 





1 1 + hx + ay 

F FORGET-ME-NOT, a Computer Code for 
the Study of Polarized Beam Dynamics 

The computer code FORGET-ME-NOT has been written for the study 
of unpolarized and polarized beam dynamics and among other things includes 
the following important options: 

1. Calculation of the strengths of the imperfection spin resonances and 
first order intrinsic spin resonances with betatron oscillations with the help 
of an averaging method. 

2. Calculation of one-turn Taylor maps for orbit and spin motion up to 
arbitrarily high order with respect to the amplitudes of the betatron and 
synchrotron oscillations and determination of 

2.1. Invariant functions of the orbit motion. 

2.2. Equilibrium polarization direction. 

2.3. Dependence of orbit and spin tunes on the invariants of orbit motion 
(spread of orbit and spin tunes). 

3. Numerical tracking of particles with spin in accelerators and storage 
rings preserving: 

3.1. Symplecticity with respect to the 6-D orbit motion. 

3.2. Orthogonality with respect to the 3-D spin motion. 

All options use the same physical model. The use of various approaches 
allows us to understand the computed results from various points of view. 
FORGET-ME-NOT, for example, has been applied to the investigation of 
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schemes for preserving the polarization in the TRIUMF KAON Booster || |7J . 
to the investigation of spin motion at high energies in the HERA proton ring 



11] and to the study of the possibility to accelerate the polarized proton 



beam in the Nuclotron ring in Dubna |32[ 
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